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ABSTRACT 

We present a study of 15 long-duration 7-ray burst (GRB) host galaxies at z > 2. The GRBs are selected with 
available early-time afterglow spectra in order to compare interstellar medium (ISM) absorption-line properties 
with stellar properties of the host galaxies. In addition to five previously studied hosts, we consider new 
detections for the host galaxies of GRB 050820 and GRB 060206 and place 2-er upper limits to the luminosities 
of the remaining unidentified hosts. We examine the nature of the host galaxy population and find that (1) 
the UV luminosity distribution of GRB host galaxies is consistent with expectations from a UV luminosity 
weighted random galaxy population with a median luminosity of (L(UV)) = 0.1 L*; (2) there exists a moderate 
correlation between UV luminosity and Sill A 1526 absorption width, which together with the observed large 
line widths of W(1526) > 1.5 A for a large fraction of the objects suggests a galactic outflow driven velocity 
field in the host galaxies; (3) there is tentative evidence for a trend of declining ISM metallicity with decreasing 
galaxy luminosity in the star-forming galaxy population at z = 2 - 4; (4) the interstellar UV radiation field is 
found w 35- 350 x higher in GRB hosts than the Galactic mean value; and (5) additional galaxies are found 
at < 2" from the GRB host in all fields with known presence of strong Mg II absorbers, but no additional 
faint galaxies are found at < 2" in fields without strong Mg II absorbers. Our study confirms that the GRB 
host galaxies (with known optical afterglows) are representative of unobscured star-forming galaxies at z > 2, 
and demonstrates that high spatial resolution images are necessary for an accurate identification of GRB host 
galaxies in the presence of strong intervening absorbers. 

Subject headings: gamma rays: bursts — ISM: abundances — ISM: kinematics — intergalactic medium 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) are among the most energetic 
events in the universe. In particular, long-duration GRBs 
are believed to originate in the catastrophic death of massive 
stars (e.g. Woosley 1993; Paczyhski 1998; Bloom et al. 2002; 
Stanek et al. 2003; Hjorth et al. 2003; see Woosley & Bloom 
2006 for a recent review). Since massive stars evolve rapidly, 
long-duration GRBs should probe instantaneous star forma- 
tion out to the highest redshifts (e.g. Wijers et al. 1998) with 
the afterglows serving as signposts to starburst galaxies in the 
distant universe. 
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Many bursts are followed by optical afterglows that can 
briefly exceed the absolute brightness of any known quasar 
by orders of magnitude (e.g. Akerlof et al. 1999; Kann 
et al. 2007; Bloom et al. 2008) and serve as bright back- 
ground sources for probing intervening gas along the line 
of sight. Early-time, high-resolution spectroscopy of GRB 
afterglows have revealed numerous absorption features pro- 
duced by ground-state and excited-state ions in the interstel- 
lar medium (ISM) of the host galaxies (e.g. Vreeswijk et al. 
2004; Prochaska, Chen, & Bloom 2006; Vreeswijk et al. 
2007; D'Elia et al. 2008). Detailed studies based on com- 
parisons of absorption-line strengths have yielded accurate 
constraints on the host ISM properties, including gas den- 
sity, temperature, chemical composition, and kinematics of 
the GRB host environment (e.g. Fynbo et al. 2006; Savaglio 
2006; Prochaska et al. 2007a). Specifically, roughly 50% of 
known GRBs at z > 2 are found in ISM of neutral gas column 
density N(U I) > 10 21 cm" 2 (Jakobsson et al. 2006; Chen et 
al. 2007a). In addition, the host ISM generally exhibit moder- 
ate chemical enrichment with a median metallicity of > 1/10 
solar, although with a substantial scatter over the range from 
1/100 to ~ 1/2 solar values (Fynbo et al. 2006a; Savaglio 
2006; Prochaska et al. 2007a). Comparisons of different 
ionic abundances also show that there exists a large differen- 
tial depletion with [Zn/Fe] > +0.6 dex and a/Fe > 0.4 dex, 
confirming the presence of a large amount of gas mass and 
a chemical enrichment history dominated by massive stars 
(Savaglio 2006; Prochaska et al. 2007a). Finally, there is a 
lack of molecular gas despite the presence of a large N(H I) 
(Fynbo et al. 2006; Tumlinson et al. 2007). 

Interpretations of these absorption-line data are not straight- 
forward. In particular, the observed low-metal content in the 
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GRB host ISM, as opposed to optically selected luminous 
starburst galaxies (e.g. Shapley et al. 2004), can be explained 
if the progenitor stars originate in the outskirts of a luminous 
galaxy or if the host galaxies are underluminous and have on 
average lower metallicities. A local-galaxy analogue of this 
is seen with the so-called luminosity-metallicity relation (e,g. 
Tremonti et al. 2004 for nearby galaxies). Recently, Fynbo 
et al. (2008) considered both scenarios and showed that the 
metallicity distribution of GRB hosts is consistent with the 
expectation that these host galaxies represent a weighted star- 
forming galaxy population according to the on-going star for- 
mation rate (SFR). In addition, the large atomic gas column 
density in contrast to the lack of molecular gas may be due 
to either an enhanced UV radiation field in the star-forming 
regions near the GRB progenitor or a relatively low dust and 
metal content indicated in the absorption-line data (e.g. Tum- 
linson et al. 2007; Whalen et al. 2008). Recent detections of 
the 2175-A dust absorption feature in GRBs 070802 (Kriih- 
ler et al. 2008; Eliasdottir et al. 2008) and 080805 (Jakobs- 
son et al. 2008) offer an important test for this scenario. Re- 
gardless, these issues have direct impact on our understand- 
ing of both the GRB progenitors and star formation physics. 
Supplemental imaging and spectroscopic observations of the 
host galaxies are necessary for accurate interpretations of the 
absorption-line measurements. 

The transient nature of optical afterglows allows deep imag- 
ing and spectroscopic studies of galaxies close to the lines 
of sight, including the hosts, when the afterglows disappear 
(c.f. M0ller et al. 2002a; Chen & Lanzetta 2003 for searches 
of damped Lya absorbing galaxies along quasar sightlines). 
At z > 2, where accurate ISM abundance measurements are 
available based on afterglow absorption-line spectroscopy, 
only four host galaxies have been unambiguously identified 
(see Savaglio et al. 2008 for a compilation). Comparison stud- 
ies between known ISM properties from afterglow absorp- 
tion spectroscopy and the observed morphology and luminos- 
ity of the host galaxies have been published individually for 
GRB 000926 (Castro et al. 2003), GRB 01 121 1 (Vreeswijk et 
al. 2006), GRB 021004 (Fynbo et al. 2005), and GRB 030323 
(Vreeswijk et al. 2004). The four GRB host galaxies together 
show an order-of-magnitude scatter in their rest-frame UV lu- 
minosity and ISM metallicity. 

At z < 2, where the majority of common ISM absorption 
features occur at ultraviolet wavelengths and spectroscopic 
observations become challenging on the ground, > 40 GRB 
host galaxies have been identified in late-time imaging follow- 
up. These known host galaxies provide important insights 
for understanding the nature of GRB progenitors at interme- 
diate redshifts. First, Chary et al. (2002) compared the es- 
timated SFR and total stellar mass for 12 GRB host galax- 
ies and found that the host galaxies have on average higher 
SFR per unit stellar mass than local starburst galaxies. Le 
Floc'h et al. (2003) examined the optical and near-infrared 
R-K colors and rest-frame B-band luminosity function of 15 
GRB host galaxies at (z) ~ 1. They found that these host 
galaxies have on average bluer colors and fainter luminosity 
((Lb) ~ 0. 1 L*) than random star-forming galaxies at the same 
redshift range. Additional mid-infrared imaging observations 
by these authors showed that this is not due to dust extinction 
(Le Floc'h et al. 2006). Similarly, Christensen et al. (2004) 
studied the optical and near-infrared spectral energy distribu- 
tions (SEDs) of ten GRB host galaxies at (z) = 0.85. They 
found based on a comparison with ~ 1000 galaxies identi- 



fied at a similar redshift range in the Hubble Deep Fields that 
GRB host galaxies have on average younger stellar age and 
shorter characteristic star-forming time scale. In addition, 
Fruchter et al. (2006) compared HST images of > 40 GRB 
host galaxies at (z) ~ 1 with the host galaxies of core-collapse 
supernovae (SNe). They found that the majority of GRB host 
galaxies exhibit irregular morphology, unlike the host galax- 
ies of core-collapse SNe. Recently, independent studies by 
Castro Ceron et al. (2008) and Savaglio et al. (2008) show 
that (z) < 1 GRB host galaxies contain on average lower stel- 
lar mass than field star-forming galaxies. Together, these re- 
sults show that GRB host galaxies at (z) ~ 1 represent a rel- 
atively young dwarf population that have experienced recent 
on-going star-forming episodes. 

While GRB host galaxies based on the intermediate- 
redshift sample appear to be underluminous and low mass 
systems, it is not clear whether the long-duration GRBs orig- 
inate preferentially in relatively metal dificient star-forming 
regions (c.f. Wolf & Podsiadlowski 2007; Modjaz et al. 2008). 
A low-metallicity environment is favored by popular progen- 
itor models so that the progenitor star can preserve high spin 
and a massive stellar core to produce a GRB (e.g. Hirschi 
et al. 2005; Yoon & Langer 2005; Woosley & Heger 2006). 
Chemical abundance measurements for z < 2 galaxies have 
been based primarily on emission line observations of associ- 
ated H II regions. A subset of the host galaxies at z < 1 have 
been followed up spectroscopically for measuring emission- 
line fluxes. A mean metallicity of roughly 1 /4 solar is found 
but with a large scatter (e.g. Sollerman et al. 2005; Modjaz 
et al. 2008; Savaglio et al. 2008). The accuracy of emission- 
line-based abundance estimates depends sensitively on the ac- 
curacy of the calibrations between different line diagnostics 
(e.g. Kewly & Dopita 2002; Skillman et al. 2003; Kennicutt 
et al. 2003). Whether or not there exists a maximum metallic- 
ity for forming a GRB remains an open question. 

We have carried out an optical and near-infrared imaging 
survey of fields around 15 GRBs at z > 2. The GRBs are se- 
lected to have early-time afterglow spectra in order to com- 
pare ISM absorption-line properties with stellar properties. 
The goal is to identify the host galaxies and constrain their 
rest-frame UV and optical luminosities. The primary objec- 
tives are (1) to quantify the luminosity distribution of the GRB 
host galaxy populations and investigate whether or not the 
GRB host galaxies trace the typical star-forming galaxies at 
high redshift and (2) to examine whether there exists a cor- 
relation between the ISM metal content and host luminosity. 
The starburst nature of GRB hosts makes this galaxy sample 
a unique laboratory for studying star formation physics and 
stellar feedback at high redshift. We adopt a ACDM cosmol- 
ogy, fiivi = 0.3 and £7a = 0.7, with a dimensionless Hubble 
constant h = // /(100 km s" 1 Mpc -1 ) throughout the paper. 

2. THE GRB SAMPLE 



We generated a sample of 15 GRBs at z G rb > 2 for studying 
the nature of high-redshift GRB host galaxies and for compar- 
ing the host properties with those of field galaxies in deep sur- 
veys. The GRBs were selected to have early-time, moderate- 
to-high resolution afterglow spectra available for measuring 
the underlying neutral hydrogen column density. In addition, 
eight of the GRBs have sufficiently high-spectral resolution 
for constraining the chemical abundances of the host ISM. 
Only five of these fields have late-time images published in 
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TABLE 1 

Summary of the Optical and Near-infrared Imaging Data 



Field RA(J2000) Dec(J2()00) z GRB log/V(HI) References 11 Instrument Filter EXPTIME (s) FWHM (") 

GRB011211 11:15:17.98 -21:56:56.2 2.140 20.4±0.2 (1) HST/STIS b Clear 19544 0.1 

Magellan/PANIC H 12960 0.4 

GRB020124 09:32:50.81 -11:31:10.6 3.198 21.7±0.2 (2) HST/STIS C Clear 24798 0.1 

Magellan/PANIC H 13860 0.4 

GRB030323 11:06:09.38 -21:46:13.3 3.372 21.90±0.07 (3) HST/ACSAVFC' 1 F606W 5928 0.1 

Magellan/PANIC H 14400 0.4 

GRB030429 12:13:07.50 -20:54:49.7 2.658 21.6±0.2 (4) Magellan/PANIC H 7740 0.4 

GRB050401 16:31:28.81 +02:11:14.2 2.899 22.6±0.3 (5) Magellan/PANIC H 6300 0.6 

Keck/LRIS g 2540 1.1 

Keck/LRIS R c 2460 1.0 

GRB050730 14:08:17.13 -03:46:16.7 3.968 22.15±0.10 (6) Magellan/MagIC 2700 0.6 

Keck/LRIS g 3900 1.2 

Keck/LRIS R c 3900 1.0 

GRB050820A 22:29:38.11 +19:33:37.1 2.615 21.0±0.1 (7) HST/ACSAVFC 6 F625W 2238 0.1 

HST/ACSAVFC F775W 4404 0.1 

HST/ACSAVFC F850LP 14280 0.1 

Magellan/PANIC H 8460 0.6 

Keck/LRIS g 2620 0.7 

Keck/LRIS R c 2500 1.9 

GRB050908 01:21:50.75 -12:57:17.2 3.343 17.55±0.10 (8) Magellan/PANIC H 12150 0.5 

GRB050922C 21:09:33.08 -08:45:30.2 2.199 21.5 ±0.1 (7) Magellan/PANIC H 12560 0.5 

Keck/LRIS g 3870 1.2 

Keck/LRIS R c 3660 1.0 

GRB060206 13:31:43.42 +35:03:03.6 4.048 20.85±0.10 (9) HST/ACSAVFC f F814W 9886 0.1 

GRB060607 21:58:50.40 -22:29:46.7 3.075 16.85±0.10 (10) Magellan/PANIC H 18540 0.5 

Keck/LRIS g 2550 1.0 

Keck/LRIS R c 2430 1.5 

GRB070721B 02:12:32.97 -02:11:40.4 3.626 21.50±0.20 (11) Magellan/PANIC H 10320 0.5 

Previously Published Fields 

GRB000301C 16:20:18.60 +29:26:36.0 2M0 21.2 ±0.5 (12) HSTTSTLP Clear 16422 0l 

GRB000926 17:04:09.62 +51:47:11.2 2.038 21.3±0.3 (13) Castro et al. (2003) 

GRB021004 00:26:54.68 +18:55:41.6 2.329 19.5±0.5 (14) Fynbo et al. (2005) 

a (1): Vreeswijk et al. (2006); (2) Hjorth et al. (2003); (3) Vreeswijk et al. (2004); (4) Jakobsson et al. (2004); (5) Watson et al. (2006); (6) Chen et al. (2005); 
(7) Prochaska et al. (2007); (8) Fynbo, private communication; (9) Fynbo et al. (2006); (10) Chen et al. (2007), see also § 4.11; (11) Malesani et al. (2007); (12) 
Jensen et al. (2001); (13) Fynbo et al. (2002); (14) Fynbo et al. (2005), but see § 4.13 for discussion. 

b The imaging data were retrieved from the HST data archive; PID = 8867. The host galaxy was clearly detected in these images that have been analyzed and 
published in Jakobsson et al. (2003). 

c The imaging data were retrieved from the HST data archive; PID = 9180. A detailed analysis of these space images is published in Berger et al. (2002). No 
emission from the host galaxy was found. 

d The imaging data were retrieved from the HST data archive; PID = 9405. The host galaxy was found in the first-epoch images (exptime = 1920 s), which have 

been analyzed and published by Vreeswijk et al. (2004). 

e The imaging data were retrieved from the HST data archive; PID = 10551. 

f The imaging data were obtained through our own program before ACS failed; PID = 10817. 

g The imaging data were retrieved from the HST data archive; PID = 8189. The host galaxy was marginally detected in these images that have been analyzed and 
published in Fruchter et al. (2006). 

the literature. To identify the stellar counterpart of the host GRBs in our sample span a range from z = 2.04 to z = 4.05 
galaxies, we have carried out a near-infrared imaging survey (left panel of Figure 1). The neutral hydrogen column densi- 
of these fields using PANIC (Martini et al. 2004) and the H fil- ties of the GRB host ISM span a range from log N(H I) = 16.9 
ter on the Magellan Baade telescope on Las Campanas, Chile. to log A^(H I) = 22.6 (right panel of Figure 1). The sample 
In addition, we have obtained and analyzed late-time deep size is restricted by the amount of available observing time, 
optical images available for some of these GRB fields from In comparison to GRBs with known redshift or A^(H I) in the 
either our own observations using the Advance Camera for literature, we show in Figure 1 that our sample is represen- 
Surveys (ACS; Ford et al. 1998) or unpublished data found tative of the spectroscopically identified GRB population in 
in the Hubble Space Telescope (HST) data archive 1 . A sum- the redshift and N(H I) parameter space. Note that three of 
mary of the fields is presented in Table 1, where we list the the 15 known GRB host absorbers do not contain high A^(H I) 
GRBs, their coordinates, redshift, and inferred neutral hydro- (log A^(H I) > 20.3) that would qualify them as a damped Lya 
gen column density, N(H I), in columns (1) through (5). The absorber (e.g. Wolfe et al. 2005). The observed low N(H I) in- 
corresponding references for the redshift and A^(H I) measure- dicates that the ISM in front of these GRBs is mostly ionized, 
ments are listed in column (6). 

t-u- r-no u . i I, • i i . , 3. IMAGING OBSERVATIONS AND DATA ANALYSIS 

This GRB host sample presented here is the largest sample 
of GRBs at zgrb > 2 for which both ISM absorption prop- To constrain the star formation and/or stellar population of 
erties and constraints on the emission properties of the host the host galaxies, we have carried out an optical and near- 
galaxies are available. It offers a unique opportunity to carry infrared imaging survey of nine GRB fields that have not been 
out a systematic study to understand the nature of starburst studied before. We have also obtained new near-infrared im- 
galaxies hosting GRBs at high redshift. The redshifts of the ages of the fields around GRB 011211 (previously studied by 

Jakobsson et al. 2003) and GRB 030323 (previously studied 

1 http://archive.stsci.edu/hst/. by Vreeswijk et al. 2004), and analyzed additional HST im- 
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FIG. 1. — Left: Redshift distribution of the 15 GRBs in our sample (shaded histogram), in comparison to the distributions of all GRBs with known redshifts 
(open histogram) and of those found by the Swift Satellite (Gehrels et al. 2004; dashed histogram). Right: Neutral hydrogen column density distribution of the 
GRBs in our sample in shaded histograms, in comparison to the distribution of known host N(H I) for GRBs found prior to July 2007 (see Chen et al. 2007a). 



ages of the field around GRB 030323 that were not included 
in Vreeswijk et al. (2004). At z < 3, near-infrared images of- 
fer valuable constraints for the intrinsic luminosity of the host 
galaxies at rest-frame optical wavelengths, while optical im- 
ages provide constraints for their rest-frame UV luminosities. 
Here we describe relevant imaging observations and data pro- 
cessing. 

3.1. Optical Images from the Hubble Space Telescope 

High spatial resolution and high sensitivity optical images 
obtained using the Wide Field Channel (WFC) of Advance 
Camera for Surveys (ACS) on board the Hubble Space Tele- 
scope (HST) are available for GRBs 030323, 050820A, and 
060206. The field surrounding GRB 030323 was observed 
under program ID 9405 (PI: Fruchter) using ACS/WFC and 
the F606W filter during July 2003 and December 2003. The 
observations were carried out in a sequence of four expo- 
sures of between 480 and 522 s each. The images were 
retrieved from the HST data archive. The field surround- 
ing GRB050820A was observed under program ID 10551 
(PI: Kulkarni) using ACS/WFC and the F625W, F775W, and 
F850LP filters during September 2005 and June 2006. The 
observations were carried out in a sequence of two to four 
exposures of between 400 and 807 s each. The images 
were retrieved from the HST data archive. The field sur- 
rounding GRB 060206 was observed under our own program 
(PID=10817; PI: Chen) using ACS/WFC and the F814W fil- 
ter during November and December 2006. The observations 
were carried out in the standard "ACS-WFC-DITHER-BOX" 
pattern of exposures between 1215 and 1256 s each. 

Individual exposures were reduced using standard pipeline 
techniques, corrected for geometric distortion using drizzle, 
registered to a common origin, filtered for deviant pixels 
based on a 5<r rejection criterion, and combined to form a 
final stacked image. A summary of the optical imaging obser- 
vations is presented in columns (7) through (10) of Table 1, 
which lists for each field the instrument and filter used, total 
exposure time, and the full width at half maximum (FWHM) 
of the median point spread function (PSF) as determined from 
point sources. We note that in addition to the GRBs listed 
in Table 1 we will include in the following analysis previ- 



ous HST imaging observations of GRB 000926 (z G rb = 2.038) 
and GRB 021004 (z G rb = 2.323) published by Fynbo et al. 
(2002), Castro et al. (2003) and Fynbo et al. (2005), respec- 
tively. 

3.2. Optical and Near-infrared Images from the Magellan 
and Keck Telescopes 

Optical images of GRB 050730 were obtained using 
MagIC on the Magellan Baade telescope in June 2008. The 
observations were carried out in one set of three exposures, 
900 s in duration. Dither offsets of 15" were applied be- 
tween exposures. The sky condition was photometric with a 
mean seeing of 0.6". Individual exposures were first corrected 
for pixel-to-pixel variation using a flat-field image formed by 
median filtering sky images obtained during even twilight. 
Fringes are not apparent in the MagIC-/' images. Next,the 
processed individual images obtained on the same night were 
registered to a common origin, filtered for deviant pixels 
based on a 5 a rejection criterion and a bad pixel mask formed 
using the flat-field frames, and stacked according to a weight- 
ing factor that is proportional to the inverse of the sky vari- 
ance. The photometric zero points were determined using five 
SDSS southern standard stars (Smith et al. 2002) observed 
during the night. 

Optical R c and g images of GRBs 050401, 050730, 
050820A, 050922C, and 060607 were obtained using the Low 
Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) on 
the the Keck I telescope in May and July of 2006. Integration 
times varied from source to source and are given in Table 1, 
but generally consisted of dithered exposures of about 600 s 
each and total integration times of 40 to 65 minutes per field. 
Individual exposures were reduced and combined using stan- 
dard techniques. Photometric calibration was performed us- 
ing a series of exposures on the Landolt field Markarian A at 
three different elevations during the July run and a single visit 
to the PG 22 1 3 field during the May run. The observation con- 
ditions were photometric. Astrometry was performed using a 
large sample of USNO B1.0 standard stars in each field. 

Near-infrared H images of GRBs 01 1211, 030323, 050401, 
050820A, and 050922c were obtained using PANIC on the 
Magellan Baade telescope in February 2004, May 2006, and 
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August 2007. The observations were carried out in five or nine 
sets of three to four exposures, 45 to 60 s in duration. Dither 
offsets of eight to 15 arc sec were applied between different 
sets of exposures in a slanted square pattern. 

Individual exposures were first corrected for pixel-to-pixel 
variation using a flat-field image formed by median filtering 
all the images obtained on the same night. Next, we corrected 
for geometric distortion in individual flat-field images using 
the IRAF geomap task, according to a 2D distortion map pro- 
vided by the PANIC instrument team. Next,the processed in- 
dividual images obtained on the same night were registered 
to a common origin, filtered for deviant pixels based on a 
5 a rejection criterion and a bad pixel mask formed using the 
flat-field frames, and stacked according to a weighting factor 
that is proportional to the inverse of the sky variance. Im- 
ages obtained during non-photometric nights were scaled to 
match the fluxes of common objects observed during photo- 
metric conditions. Individual stacked images from different 
nights were combined to form a final image of each field. The 
photometric zero points were determined using three to five 
IR standard stars (Persson et al. 1998) observed under photo- 
metric conditions. A summary of the near-infrared imaging 
observations is presented in Table 1 . 

3.3. Astrometry 

An accurate astrometric solution is necessary for the im- 
ages obtained at late times, in order to correctly identify the 
host galaxies of the GRBs. To obtain an accurate astrometric 
solution for each final stacked image, we first calibrate the as- 
trometry using < 12 USNO stars with a low-order polynomial 
fit. Next, we refine the astrometric solution using > 2 2MASS 
stars in the image by adjusting field offsets and rotation. We 
find that the final astrometric solution is accurate to am r.m.s. 
of 0.2". 

4. DESCRIPTION OF INDIVIDUAL FIELDS 

Including previously published images around GRBs 
00030 1C, 000926 and 021004, we have now collected deep 
optical and/or near-infrared images surrounding 15 GRBs at 
z > 2. All of these GRBs have early-time afterglow spec- 
tra available that allow us to determine both the gas proper- 
ties of the GRB host galaxies and the line-of-sight proper- 
ties of absorbers foreground to the hosts. We will show that 
known line-of-sight properties from afterglow spectra are im- 
portant for accurately identifying the GRB host galaxies (e.g. 
GRB 030429, GRB 060206, GRB070721B, and see Pollack 
et al. 2008 for a rich galaxy field around GRB 060418). Here 
we describe contraints on the emission properties of individ- 
ual host galaxies from available imaging data, together with 
a brief description of known absorption-line properties from 
afterglow spectroscopy. A summary of the emission and ab- 
sorption properties of these z > 2 GRB host galaxies is pre- 
sented in Table 3. 

4.1. GRB011211 atz = 2.U2 

This burst was detected by BeppoSAX and the optical tran- 
sient was found w 10 hours later with R = 19 (Grav et al. 
2001). Low-resolution (FWHM = Sv rj 680 km s" 1 ) optical 
spectra of the afterglow were taken by different groups, which 
revealed multiple absorption features indicating a source red- 
shift of z = 2.142 (Fruchter et al. 2001; Gladders et al. 2001). 
Analyzing existing afterglow spectra obtained using FORS2 
on the VLT telescopes, Vreeswijk et al. (2006) reported a total 
neutral hydrogen column density of log A^(H I) = 20.4 ± 0.2. 
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Fig. 2. — Registered optical (left three panels) and near-infrared H (right 
panel) images of the field around GRB 01 121 1 at zqrb = 2.142 after the after- 
glow had faded. The epoch during which the images were taken is indicated 
at the bottom of each panel. The host appears to consist of three compact 
regions, all of which are confirmed to be at the host redshift based on the 
presence of Lya emission in ground-based narrow-band imaging follow-up 
(Fynbo et al. 2003a). The GRB is found by Jakobsson et al. (2003) to orig- 
inate in the fainter region, southwest of the center blob, as indicated by the 
fading optical transient (arrow in the left panel). The mean FWHM of the 
PSF in the H image was found to be f» 0.37" based on an average of 10 stars 
across the PANIC field. Two of the compact regions are detected in the H- 
band image, but the region directly associated with the burst does not exhibit 
observable flux. We measure an //-band magnitude of AB(H) = 25.0 ±0.3 
over a 1" diameter aperture for the host galaxy. 

In addition, the authors applied a curve-of-growth (COG) 
analysis over a series of absorption features found in the af- 
terglow spectra and derived [Si/H] = -0.9^ 4 and [Fe/H] = 
-1 .3 ± 0.3 2 . The large A^(H I) indicates that even the weakest 
absorption features detected in the afterglow spectrum may 
be saturated due to the low resolution and low S/N of the 
data. Here we consider the abundance measurements based 
on the COG analysis lower limits to the intrisic ISM metallic- 
ity and adopt [Si/H] > -1.3 for the ISM of the host galaxy. 
However, we note that the absence of relatively weak tran- 
sitions such as Sill A 1808 in the afterglow spectrum indi- 
cates [Si/H] < -0.7. Finally, this line of sight exhibits no 
Mgll absorbers of rest-frame absorption equivalent width 
W(2796) > 1 A between z = 0.359 and z = 2 (Prochter et al. 
2006). 

Imaging follow-up of the field around GRB 01 1211 was 
carried out during four different epochs with STIS and the 
clear filter on board HST The imaging data were reduced 
and analyzed by Jakobsson et al. (2003), who reported the 
host galaxy has an /?-band magnitude of AB(R) = 25.15 ±0.11 
over a 1 "-diameter aperture. In addition, Lya emission 
was detected in a ground-based narrow-band imaging sur- 
vey by Fynbo et al. (2003a), who measured a total flux of 
/(Lya) = (2.8 ±0.8) x 10" 17 erg s _I cm -2 and derived an SFR 
of 0.8 ±0.2 Mq yr" 1 . 

We have observed this field using PANIC on Magellan in 
February 2004, and obtained a total integration of 216 min- 
utes. The mean FWHM of the PSF was found to be ss 0.37" 
based on an average of 10 stars across the PANIC field. The 
final stacked image is presented in Figure 2, together with 
HST STIS images obtained roughly 14, 26, and 32 days af- 
ter the initial burst. The optical and H images are registered 
to a common origin. The host galaxy is marked by a circle 

2 Chemical abundances are measured relative to solar values and defined 
as [M/H] = log(M/H)-log(M/H) . 
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FIG. 3. — Schematic diagram to illustrate the photometric zero point offset 
neccessary to be included for sources at z = 3.2, from a clear bandpass that 
covers a spectral range over A = 2000- 10, 000 A to a typical Johnson R filter. 
The historgram represents a typical starburst spectral template with additional 
absorption at A res t < 1215 A due to the IGM Lya forest and the ISM internal 
to the host. 

of 1" radius. The host appears to consist of three compact 
(presumably star-forming) regions in the STIS images, all of 
which are confirmed to be at the host redshift based on the 
presence of Lya emission in the ground-based narrow-band 
images presented in Fynbo et al. (2003a). The GRB is found 
to originate in the faintest of the three blobs, southeast of the 
center one. Two of the compact regions are detected in the H- 
band image, but the region directly associated with the burst 
does not exhibit any detectable flux. We measure an //-band 
magnitude of AB(H) = 25.0 ± 0.3 over a 1" diameter aperture 
for the host galaxy. 

At z = 2.142, the observed /?-band magnitude corresponds 
to a rest-frame absolute magnitude at 2000 A of M(2000) - 
5 log h = -19.2 ±0.1. The observed //-band magnitude al- 
lows us to derive a rest-frame absolute B-band magnitude of 
M AB {B)-5 log h = -19.1 ±0.3 for the GRB host galaxy. 

4.2. GRB 020124 at zgrb = 3.198 

This burst was detected by HETE (Ricker et al. 2002) and 
the optical transient was found rj 22 hours later with R ~ 18.5 
(Price et al. 2002). Low-resolution (5v ks 680 km s" 1 ) op- 
tical spectra of the afterglow were obtained using FORS1 
on the VLT Melipal telescope, which revealed multiple ab- 
sorption features indicating a source redshift of z = 3.198 
(Hjorth et al. 2003). The total neutral hydrogen column den- 
sity derived based on the observed Lya absorption line is 
log JV(H I) = 21.7 ± 0.2 (Hjorth et al. 2003). The low reso- 
lution and low S/N (« 3-4) of the data did not allow an ac- 
curate measurement of the chemical content in the host ISM. 
No information is available for the line-of-sight properties of 
additional intervening absorbers. 

Imaging follow-up of the field around GRB 020124 was 
carried out roughly 18 and 21 days after the burst with STIS 
and the clear filter on board HST The imaging data were re- 
duced and analyzed by Berger et al. (2002). While the OT 
was detected in the first epoch image, no detectable flux was 
found at the OT position in the second epoch image. Berger 
et al. (2002) placed an upper limit for the /?-band magnitude 
of the host galaxy at R > 29.5. We have retrieved the imaging 
data from the HST archive and determined a 5-er limit of the 
second epoch image at AB(clear) = 29.4 over a 0.5" diame- 




FlG. 4. — Registered optical (left two panels) and near-infrared H (right 
panel) images of the field around GRB 020124 at zgrb = 3.198. The epoch 
during which the images were taken is indicated at the bottom of each panel. 
The //-band image has been smoothed using a Gaussian kernel of FWHM = 
0.4", which is roughly the size of the PSF. While the OT is still visible in 
the first-epoch image (left panel), the host is not detected in either of the two 
late-time images. We place a 2-cr //-band limiting magnitude of AB(H) = 26. 1 
over a 0.5" diameter aperture for the host galaxy. 

ter aperture. To derive the corresponding detection limit, we 
take into account the bandpass difference between STIS/clear 
and R and the IGM opacity that reduces most of the light at 
Aobs < 4000 A for sources at z = 3.198 (Figure 3). We apply 
a 0.4 mag offset in the photometric zero point and conclude 
that the host galaxy is fainter than AB(R) = 29 at the 5-er level 
of significance. 

We have also observed this field using PANIC on Mag- 
ellan in February 2004, and obtained a total integration of 
231 minutes. The mean FWHM of the PSF was found to be 
k 0.4". The final stacked image is presented in Figure 4, to- 
gether with available HST STIS images. The optical and H 
images are registered to a common origin. The //-band image 
has been smoothed using a Gaussian kernel of FWHM = 0.4". 
The position of the GRB is marked by a circle of 0.5"radius. 
No detectable flux is found at the location of the GRB, but 
some emission features are observed with AB(H) w 26.4 at 
A 9 = 0.4" angular distance away. We place a 2-cr limit in the 
observed //-band magnitude of AB(H) = 26.1 over a 0.5" di- 
ameter aperture for the host galaxy. At z = 3.198, the observed 
//-band magnitude limit allows us to derive a limiting rest- 
frame absolute magnitude of Mab(3960)-5 log h > -18.8 for 
the GRB host galaxy. 

4.3. GRB 030323 at z G rb = 3.372 

This burst was detected by HETE-2 (Graziani et al. 2003) 
and the optical transient was found « 7.6 hours later with R = 
18.7 (Gilmore et al. 2003a). Moderate-resolution (Sv w 150 
km s" 1 ) optical spectra of the afterglow, covering A = 4560- 
73 10 A, were taken using FORS2 on the VLT Yepun telescope 
by Vreeswijk et al. (2004), revealing multiple absorption fea- 
tures from ions in both ground states and excited states that 
are consistent with a source redshift of z = 3.372. These au- 
thors estimated the total neutral hydrogen column density of 
log JV(H I) = 2 1.90 ±0.07 in the host ISM and chemical abu- 
dances of [S/H] =-1.26±0.2 and [Fe/H] =-1.47 ±0.11. 
The observed large column densities of various ions sug- 
gest that these reported values represent only lower limits to 
the intrinsic abundances of these ions. We therefore adopt 
[S/H] > -1.26 for the ISM of the host galaxy. This line of 
sight exhibits no strong Mgll absorbers at z = 0.824- 1.646 
(Prochter et al. 2006). 

Imaging follow-up of the field around GRB 030323 was 
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FIG. 5. — Registered optical (left panel) and near-infrared H (right panel) 
images of the field around GRB 030323 at zgrb = 3.372. The epoch during 
which the images were taken is indicated at the bottom of each panel. The 
host galaxy is identified at AS = 0.22" from the position of the OT with 
AB(F606W) = 27.4 ± 0.1 over a 0.5" diameter aperture. The //-band im- 
age has been smoothed using a Gaussian kernel of FWHM = 0.4" , which 
is roughly the size of the PSF. The host is not detected in the H image. We 
place a 2-cr //-band limiting magnitude of AB(H) = 26.2 over a 0.5" diameter 
aperture for the host galaxy. 

carried out during two different epochs with ACS and the 
F606W filter on board HST. The first epoch imaging data 
with a total exposure time of 1920 s were reduced and an- 
alyzed by Vreeswijk et al. (2004), who identified the host 
galaxy at A (9 = 0.14" from the position of the OT and mea- 
sured AB(F606W) = 28.0 ± 0.3 over a 0.3" diameter aperture. 
Additional imaging data were obtained five month later with 
HST using the same instrument setup. A stack of all avail- 
able imaging data from the HST data archive shows a clear 
detection of the host galaxy at A9 = 0.22" from the posi- 
tion of the OT. This position is consistent with the position 
of Vreeswijk et al. to within the astrometric uncertainties. Af- 
ter correcting for the Galactic extinction (E(B-V) = 0.049 
according to Schlegel et al. 1998), we measure a total flux of 
AB(F606W) = 27.4 ± 0. 1 over a 0.5" diameter aperture. 

We have also observed this field using PANIC on Magellan 
in February 2004, and obtained a total integration of 240 min- 
utes. The mean FWHM of the PSF was found to be w 0.4". 
The final stacked image is presented in Figure 5, together with 
a stack of available ACS/F606W images. The optical and H 
images are registered to a common origin. The host galaxy 
is marked by a circle of 0.5" radius. The host appears to be 
extended in the ACS image, but not detected in the stacked 
H image. We place a 2-cr limit in the observed //-band mag- 
nitude of AB(H) = 26.2 over a 0.5" diameter aperture for the 
host galaxy. 

At z = 3.372, the observed F606W magnitude corresponds 
toM /tB (1400)-5 log h = -17.6 ±0.1 for the GRB host galaxy. 
The observed //-band magnitude limit allows us to derive 
a limiting rest-frame absolute magnitude of Mab(3800) - 
5 log h > -18.8 for the GRB host galaxy. 

4.4. GRB 030429 at z G rb = 2.658 

This burst was detected by HETE-2 (Doty et al. 2003) and 
the optical transient was found « 3.5 hours later (Gilmore 
et al. 2003b) at A 6 « 1 .2" southeast of an extended source 
of R w 24 (Fynbo et al. 2003b). Low-resolution (Sv « 680 
kms" 1 ) optical spectra of the afterglow and the extended 
source were obtained using FORS 1 on the VLT Melipal tele- 
scope. The spectrum of the afterglow exhibits multiple ab- 
sorption features that are consistent with z = 2.658 (Jakobsson 
et al. 2004). In contrast, the spectrum of the extended source 




FIG. 6. — Near-infrared H image of the field around GRB 030349 at z G RB = 
2.658. The mean FWHM of the PSF is 0.35". A smooth version of the H- 
band image is presented in the right panel. No detectable flux is found at the 
immediate location of the GRB, but extended emission features (pointed by 
an arrow in the left panel) are clearly visible at as 1 " southwest of the OT, or 
1.3" south of the galaxy at z = 0.841. We measure AB(H) = 20.57 ± 0.05 for 
the foreground galaxy, and AB(H) = 24.4 ± 0. 1 for the host candidate. 

exhibits a single emission feature at A = 6858 A. Jakobsson 
et al. (2004) identified this emission as [Oil] at z = 0.841, 
which is supported by the presence of an absorption feature 
in the afterglow spectrum if interpreted as a Mg II absorber at 
z = 0.841. The GRB host galaxy therefore remains unidenti- 
fied. These authors also reported a total neutral hydrogen col- 
umn density of log A^(H I) = 21.6 ±0.2 based on the observed 
Lya absorption feature. The low resolution of the spectrum 
did not allow an accurate measurement of the chemical con- 
tent in the host ISM. 

We have observed the field around GRB 030429 using 
PANIC on Magellan in February 2004, and obtained a total 
integration of 129 minutes. The mean FWHM of the PSF was 
found to be m 0.35". The final stacked image is presented in 
Figure 6, together with a smoothed version with a Gaussian 
kernel of FWHM = 0.4". The position of the GRB is marked 
by a circle of 0.5"radius. No detectable flux is found at the 
immediate location of the GRB, but additional emission fea- 
tures are present at ss 1" southwest of the OT, or 1.3" south 
of the galaxy at z = 0.841 . We measure AB(H) = 20.57 ± 0.05 
for the foreground galaxy, and AB(H) = 24.4 ± 0. 1 for the host 
candidate. 

The emission morphology resembles those found for the 
hosts of GRB 000926 and GRB 011211, with the OT originat- 
ing in the faintest of the three emission blobs. The close prox- 
imity to the foreground galaxy suggests that the host may be 
gravitationally lensed by this foreground object, although the 
large angular separation implies that the foreground galaxy 
would have to be rj 4 x as massive as the Milky Way (see also 
Jakobsson et al. 2004). Additional images at optical wave- 
lengths are necessary to investigate this lensing hypothesis 
further. 

If the extended feature is associated with the host galaxy, 
then the OT would be at a projected distance of p = 5.4 h 
kpc from the center of the galaxy. In the following discus- 
sion, we consider the observed brightness as an upper limit 
to the brightness of the GRB host and derive a rest-frame B- 
band magnitude of M AB (B)-5 log h > -20. 1 for the GRB host 
galaxy. 

4.5. GRB 050401 at z G rb = 2.899 

This burst was detected by Swift (Barbier et al. 2005) and 
prompt localization of the source was reported by Angelini 
et al. (2005) from the X-ray afterglow. The optical transient 
was found w 1 hour later with R = 20.3 (McNaught & Price 



8 



2005). Low-resolution (Sv w 680 km s _I ) optical spectra of 
the afterglow were taken using FORS2 on the VLT Antu tele- 
scope (Watson et al. 2006), revealing multiple absorption fea- 
tures from ions in both ground states and excited states that 
are consistent with a source redshift of z = 2.899. Watson et 
al. (2006) estimated the total neutral hydrogen column den- 
sity of logJV(H I) = 22.6 ±0.3 in the host ISM and metal- 
licity of [Zn/H] = -1.0 ±0.4. Similiar to the other GRB 
hosts where only low-resolution afterglow spectra are avail- 
able, the observed large column densities of various ions sug- 
gest that these reported values represent only lower limits to 
the intrinsic abundances of these ions. We therefore adopt 
[Zn/H] > -1.0 for the ISM of the host galaxy. An additional 
set of ionic transitions (such as Al II and Fe II) was also re- 
ported at z = 2.5 by Watson et al. (2006). 

We have observed the field around GRB 050401 using 
PANIC on Magellan in May 2006, and obtained a total in- 
tegration of 105 minutes. The mean FWHM of the PSF was 
found to be ps 0.6". No detectable flux is seen at the immedi- 
ate location of the GRB. We place a 2-a limit in the observed 
H-band magnitude of AB(H) = 25. 1 over a 1" diameter aper- 
ture for the host galaxy. At z = 2.899, the observed //-band 
magnitude limit allows us to derive a limiting rest-frame ab- 
solute magnitude of M AB (B)-5 log h > -19.6 for the GRB 
host galaxy. 

Optical images of the field around GRB 050401 have also 
been obtained using LRIS and the g and R c filters on the Keck 
I telescope in May 2006. The mean FWHM of the PSF was 
found to be w 1.1" in the combined g image and w 1.0" in 
the combined R c image. The stacked images are presented in 
Figure 7, which have been smoothed with a Gaussian kernel 
of FWHM =1". The position of the GRB is marked by a cir- 
cle of 2" radius. At the location of the afterglow reported by 
(McNaught & Price 2005), we detect faint emission in both 
images at roughly 2-a significance level. We estimate a total 
brightness of AB(g) = 27.5 ±0.5 and AB(R C ) = 27.3 ±0.4 over 
a 2" diameter aperture for the object. Due to the presence of 
a strong absorber at z = 2.5, the identification of the observed 
faint emission is uncertain. We note, however, that at z = 2.899 
the IGM Lya forest is expected to absorb a large fraction of 
flux in the g band, resulting in g-R c ~ 0.2-0.4 mag for a flat 
spectrum source. Additional HST imaging data are necessary 
to confirm the host identification. At z = 2.899, the observed 
/? c -band magnitude would imply a rest-frame absolute mag- 
nitude of M AB (1600)-5 log h = -11. A ±0.4 for the candidate 
host galaxy. 

4.6. GRB 050730 at z G rb = 3.968 

This burst was detected by Swift (Holland et al. 2005). An 
optical transient was found promptly using the Ultraviolet- 
Optical Telescope (UVOT) on board Swift with V = 17.6 
about three minutes after the burst trigger (Holland et al. 
2005). We obtained an echelle spectrum of the afterglow, 
using the MIKE spectrograph (Bernstein et al. 2003) on the 
Magellan Clay telescope, four hours after the initial trig- 
ger. Descriptions of the data were presented in Chen et 
al. (2005) and Prochaska et al. (2007b). The spectrum 
covers a full spectral range from 3300 A through 9400 A 
with a spectral resolution of 5 v w 10 km s -1 at wavelength 
A = 4500 A and 5 v w 12 km s" 1 at A = 8000 A. The host 
of the GRB exhibits a strong damped Lya absorption fea- 
ture with log N(H I) = 22.15 ± 0.05, and metallicity [S/H] = 
-2.26 ±0.1 and [S/Fe] > +0.24 ±0.11 (see also Starling et 




FIG. 7. — Optical g (left panel) and R c (right panel) images of the field 
around GRB 050401 at z G rb = 2.899. The mean seeing is 1.1" and 1.0" in g 
and R c , respectively. At the location of the afterglow reported by (McNaught 
& Price 2005), we detect faint emission in both images at roughly 2-a sig- 
nificance level. The circle indicates a 2" angular radius around the afterglow 
position. 

al. 2005; D'Elia et al. 2007). The low metallicity and low 
a-element enhancement implies a nearly dust free medium. 
No trace of H2 is found despite the large A^(H I). The high- 
resolution, high S/N echelle data allow us to place a sensi- 
tive limit on the molecular fraction of the host ISM at /h, = 

2N(U2)/[N(UI) + 2N(H2)] < 10" 71 (Tumlinson et al. 2007). 
Additional strong absorbers are found at z = 3.56, z = 3.02, 
z = 2.25, and z = 1 .77 (Chen et al. 2005). 

We have observed the field around GRB 070530 using 
MagIC and the i' filter on Magellan in June 2008, and ob- 
tained a total integration of 45 minutes. The mean FWHM 
of the PSF was found to be w 0.6". The final stacked image 
is presented in Figure 8, which has been smoothed using a 
Gaussian kernel of FWHM = 0.6". The position of the GRB 
is marked by a circle of 1" radius. No detectable flux is seen 
at the immediate location of the GRB. We place a 2-a limit 
in the observed /'-band magnitude of AB(i') = 26.6 over a 1" 
diameter aperture for the host galaxy. At z = 3.968, observed 
/'-band magnitude limit allows us to derive a limiting rest- 
frame absolute magnitude of M4b(1500)-5 log h > -18.6 for 
the GRB host galaxy. 

Optical images of the field around GRB 050730 have also 
been obtained using LRIS and the g and R c filters on the Keck 
I telescope in May 2006. The mean FWHM of the PSF was 
found to be ps 1 .2" in the combined g image and « 1 .0" in the 
combined R c image. No detectable flux is seen at the immedi- 
ate location of the GRB. We place 2-a limits of AB(g) = 26.6 
and AB(R C ) = 26.4 over a 2" diameter aperture for the host 
galaxy. Taking into account that at z = 3.968 the IGM Lya 
forest and the large amount of neutral gas in the host ISM ab- 
sorb a large fraction of flux in the g and R c bands, we find that 
the flux limits derived from the LRIS images are consistent 
with the flux limit seen in the MagIC /'-band data. 

4.7. GRB 050820A at zgrb = 2.615 

This burst was detected by Swift (Page et al. 2005a). An 
optical transient, reported less than 1 hour after the GRB, was 
identified in data taken shortly after the trigger (Fox & Cenko 
2005; Vestrand et al. 2006). High-resolution (Sv « 7- 10 
km s" 1 ) echelle spectra of the afterglow obtained shortly af- 
ter the burst are available from both our own observations us- 
ing HIRES (Vogt et al. 1994) on the Keck I telescope and 
the ESO data archive for UVES (D'Odorico et al. 2000). 
These two spectra together cover a full spectral range from 
3300 A through 10,000 A, allowing accurate estimates of 
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FIG. 8.— A smoothed i' image of the field around GRB 050730 at z G RB = 
3.968. No detectable flux is seen at the location of the OT. We place a 2-cr 
/'-band limiting magnitude of AB(H) = 26.6 over a 1" diameter aperture for 
the host galaxy. 

chemical abudances in the ISM of the GRB host galaxy. 
Based on multiple absorption features from both ground- 
state and excited-state ions, we determine a source redshift 
of z = 2.6147 and a total neutral hydrogen column density 
of logNCH I) = 21.0 ±0.1 (Prochaska et al. 2007b; Ledoux 
et al. 2005). An absorption-line analysis of various ions 
shows that [S /H] = -0.63 ± 0.1 1 and [S /Fe] = +0.97 ± 0.09 
(Prochaska et al. 2007b), from which we derive a dust-to-gas 
ratio that is comparable to what is seen in the Small Mag- 
ellanic Cloud (SMC). Adopting the SMC dust-to-gas ratio 
(Gordon et al. 2003), we estimate a visual extinction in the 
host ISM of A v w 0.08. Finally, no trace of H 2 is found de- 
spite the large A^(H I) and moderate metallicity. The high- 
resolution, high S/N echelle data allow us to place a sen- 
sitive limit on the molecular fraction of the host ISM at 
/h 2 = 2N(H2)/[jV(HI) + 2jV(H 2 )] < 10" 6 - 5 (Tumlinson et al. 
2007). 

Imaging follow-up of the field around GRB 050820A was 
carried out with ACS and the F625W, F775W, F850LP fil- 
ters on board HST in two epochs, roughly 37 days, and nine 
months after the burst (PID=1055 1). We retrieved the imaging 
data from the HST data archive and processed the images fol- 
lowing the descriptions in § 3. Stacked images from the two 
epochs are presented in Figure 9. The OT was clearly detected 
in the first epoch image, but faded in the second epoch im- 
age which reveals faint extended emission features of the host 
galaxy. In addition to the extended low-surface feature seen at 
the position of the OT, which is identified as the host galaxy, 
we identify two compact sources at A 9 = 1 .3" north of the OT 
(Object A) and 0.4" south of the OT (Object B) from the after- 
glow lines of sight. We measure isophotal magnitudes of the 
host and Objects A and B. The isophotal apertures are defined 
based on the extent of objects found in a "white light" image, 
which is a stack of the F625W, F775W, and F850LP images. 
This "white light" image has the optimal S/N for recovering 
faint emission features and allows us to determine a common 
aperture for every object across different bandpasses. We cor- 
rect the observed brightness to account for intrinsic absorption 
in the Milky Way (E(B-V) = 0.044 along the line of sight 
according to Schlegel et al. 1998). The photometric measure- 
ments are presented in Table 2. 

We have also observed this field using PANIC on Magel- 
lan in August 2007, and obtained a total integration of 141 
minutes. The mean FWHM of the PSF was found to be 



w 0.5". The final stacked image is presented in Figure 10, 
together with the second-epoch ACS/F850LP image. The op- 
tical and H images are registered to a common origin, and 
the //-band image has been smoothed using a Gaussian ker- 
nel of FWHM = 0.5". The position of the GRB is marked 
by a circle of 0.5" radius. Extended emission at the posi- 
tion of the host galaxy is detected in the H image at » 3 a 
level of significance. We measure an //-band magnitude of 
AB(H) = 25.3 ± 0.3 within the isophotal aperture defined for 
the host galaxy in the "white light" image described above. 

At z = 2.6147, the observed F775W magnitude of the host 
galaxy corresponds to a rest-frame absolute magnitude at 
2000 A of M(2000) - 5 log h = -18.2 ± 0.06. The observed 
//-band magnitude allows us to derive a rest-frame absolute 
B-band magnitude of M AB (B)-5 log h = -19.2 ± 0.3 for the 
GRB host galaxy. Adopting Ay = 0.08 estimated above from 
absorption-line abundance ratios and the SMC extinction law 
(Gordon et al. 2003), we derive extinction corrected rest- 
frame absolute magnitudes of Mab{B)-5 log h = -19.3 ±0.3 
and M(2000)-5 log h = -18.5 ±0.06. 

Neither Object A or B exhibits detectable flux in the H 
image. We place a 2-cr limit of AB(H) > 26 for the two 
sources. The //-band photometric measurements are also pre- 
sented in Table 2. The observed optical and near-infrared 
colors of Objects A and B are relatively bluer than those of 
the GRB host galaxy and inconsistent with the expectations 
for z > 2 star-forming objects. This suggests that these are 
likely foreground galaxies. At these small angular separations 
(AO < 1.4), objects at z < 2.6 have projected distances of 
p < 8 h~ l kpc to the afterglow line of sight, and are expected 
to imprint strong Mgll A A 2796, 2803 absorption features in 
the afterglow spectrum (e.g. Chen & Tinker 2008). 

Incidentally, two strong Mg II absorbers are found at z = 
0.692 and z = 1.430 o with W(2796) = 2.99 ± 0.03 A and 
W(2796) = 1.9 ±0.1 A in the rest frame, respectively. The 
z = 0.692 Mg II absorber exhibits a complex kinematic profile 
with multiple absorption components spreading over a line- 
of-sight velocity interval of Av w 500 kms" 1 , and a non- 
negligible amount of Ca + ions (see Figure 13 in Prochaska et 
al. 2007b). The complex kinematic profile of the absorber and 
the presence of Ca + together may be explained by a sightline 
passing through an interacting system, simliar to the Milky 
Way and the Magellanic Stream (e.g. Gibson et al. 2000; Put- 
man et al. 2003). Attributing both Objects A and B to the 
z = 0.692 absorber implies a total absolute fi-band magnitude 
of M AB (B)-5 log h = -16.27 ± 0.05. Follow-up near-infrared 
Ha spectroscopy is necessary to provide conclusive identifi- 
cations of these two sources. 

Optical images of the field around GRB 050820A have also 
been obtained using LRIS and the g and R c filters on the Keck 
I telescope in July 2006. The mean FWHM of the PSF was 
found to be « 0.7" in the combined g image and w 1 .9" in the 
combined R c image. The effective seeing was significantly 
compromised in the R c images due to a focus problem with 
the LRIS red-side. Faint emission is clearly detected at the 
position of the afterglow, but we are unable to obtain accu- 
rate measurements of the host magnitude due to contaminat- 
ing light from galaxies A and B (Figure 9) in the ground-based 
images. 

For the host galaxy, the optical and //-band photometric 
measurements presented in Table 2 allow us to examine the 
star formation history and constrain the stellar mass based on 
the observed spectral energy distribution (SED). To constrain 
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Fig. 9. — Images of the field around GRB 050820A at zgrb = 2.6147. The first-epoch image (left panel), taken ss 37 days after the burst with ACS and the 
F850LP filter, clearly shows the OT (indicated by the arrow). The second-epoch image (right panel), taken nine months later, reveals extended faint blue emission 
of the host galaxy. The false-color image was formed by combining stacks of F625W, F775 W, and F850LP images. The yellow cross and circle mark the position 
of the x-ray afterglow and associated error reported by Page et al. (2005b). In addition to the host galaxy, we also point out two compact sources A and B at 
AO = 1.3" and 0.4", respectively, from the afterglow lines of sight. 
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FIG. 10.— Registered ACS F850LP image (left panel) and near-infrared H (right panel) image of the field around GRB 050820A at z G RB = 2.6147. The 
//-band image has been smoothed using a Gaussian kernel of FWHM = 0.5", which is roughly the size of the PSF. The host is detected in the H image with 
AB(H) = 25.3 ± 0.3. Neither Object A or B exhibits faint emission in the H image. We place a 2-cr limit of AB(H) > 26 for the two sources. 

of the galaxy is shown in the top panel together with the best- 
fit model. It is clear that the //-band photometry provides 
the necessary measurement for constraining the stellar popu- 
lation based on the 4000- A flux decrement. The bottom panel 
of Figure 10 shows the likelihood distribution function ver- 
sus stellar age, indicating that the last major episode of star 
formation occurred at w 400-700 Myr ago. 

Adopting the best-fit stellar synthetic model, we infer a total 
stellar mass of M* = 0.9^ 4 x 10 9 h~ 2 M . We argue that our 
inferred stellar mass is accurate, despite an absence of rest- 



the underlying stellar population, we consider a suite of syn- 
thetic stellar population models generated using the Bruzual 
& Chariot (2003) spectral library. We adopt a Salpeter initial 
mass function with a range of metallicity from 1/5 solar to 2 x 
solar and a range of star formation history from a single burst 
to exponentially declining SFR of e-folding time r = 300 Myr 
or 1 Gyr. We include intrinsic dust extinction that follows 
the SMC extinction law. Comparing the observed SED with 
model predictions allows us to constrain the stellar age. The 
results are presented in Figure 11, where the observed SED 
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FIG. 11. — Constraints on the galaxy age and star formation history. 
Top: The broad-band spectral energy distribution of the host galaxy of 
GRB050820A at z = 2.6147 (points with errorbars), in comparison to the 
best-fit template (solid histogram) and the predicted broad-band photometric 
points (open squares). The best-fit SED is charaterized by a declining SFR 
of e-folding time r = 300 Myr, 1/5 solar metallicity, and intrinsic dust ex- 
tinction of E(B-V) = 0.05 following the SMC extinction law. Bottom: The 
likelihood function of the stellar age for the host galaxy, suggesting that a 
recent starburst occurred about 400 — 700 Myr prior to the GRB explosion. 

frame near-infrared flux measurements. This is supported by 
previous studies, which show that stellar masses determined 
based on the observed 4000-A flux decrement are consistent 
with those determined based on rest-frame near-infrared lu- 
minosity to within the uncertainties (cf. Chen & Marzke 2004 
and Yan et al. 2004 for red galaxies identified in the Hubble 
Ultra Deep Field and see also Shapley et al. 2005 for a detailed 
comparison based on 72 star-forming galaxies at (z) = 2.3). 
The inferred stellar mass for the host galaxy of GRB 050820A 
at z = 2.6147 is comparable to those derived for z ~ 1 GRB 
host galaxies (e.g. Castro Ceron et al. 2006, 2008; Savaglio 
et al. 2008) but falls in the bottom 2% of the UV luminous 
galaxies studied by Shapley et al. (2005). 



4.8. GRB 050908 at z GRB = 3.344 

This burst was detected by Swift (Goad et al. 2005) and 
the optical afterglow was reported by Torii (2005) 14 minutes 
after the trigger to have R « 18.8. We obtained moderate- 
resolution (S v « 40- 150 km s" 1 ) optical spectra of the after- 
glow, using GMOS on the Gemini north telescope (Foley et 
al. 2005) and DEIMOS on the Keck II telescope (Prochaska 
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FIG. 12. — A smoothed H image of the field around GRB 050908 at z G RB = 
3.344. No detectable flux is seen at the location of the OT. We place a 2-cr 
//-band limiting magnitude of AB(H) = 26 over a 0.5" diameter aperture for 
the host galaxy. 

et al. 2005). The two spectra together provides contiguous 
spectral coverage over A = 5070-9000 A (Chen et al. 2007b). 
The afterglow spectra exhibit a suite of absorption features, 
consistent with a source redshift of z = 3.3437. Contrary 
to the majority of GRB host galaxies, the line-of-sight to- 
ward GRB 050908 displays only a modest amount of neu- 
tral gas in the host galaxy. Our initial estimate based on 
the observed Lya absorption line suggest log A^(H I) ~ 19.2. 
However, an additional spectrum obtained by Fugazza et al. 
(2005) and analyzed by Fynbo and collaborators displays 
non-negligible flux at A b s < 4000 A (corresponding to rest- 
frame wavelength range A rest < 912 A at z = 3.34), indicating 
logjV(HI) = 17.55 ±0.1 (Johan Fynbo 2007, private com- 
munication). This source repressents one of a few GRB 
sightlines observed so far that do not pass through neutral 
gas clouds in the host galaxies. Other known sources are 
GRB 021004 at z GRB = 2.329 with logJV(H I) = 19.5 ± 0.5, 
GRB 060526 at z GRB = 3.221 with log N(U I) = 20.00 ± 0.15 
(Jakobsson et al. 2006), and GRB 060607 at z GRB = 3.075 with 
logJV(H I) = 16.85±0.10 (see further discussion in §§ 4.11 
and 4. 13). A strong Mg II absorber is found at z = 1 .548 in the 
afterglow spectrum (Prochter et al. 2006). 

We have observed the field around GRB 050908 using 
PANIC on Magellan in August 2007, and obtained a total in- 
tegration of 203 minutes. The mean FWHM of the PSF was 
found to be w 0.5". The final stacked image is presented in 
Figure 12, which has been smoothed using a Gaussian kernel 
of FWHM = 0.5". The position of the GRB is marked by a cir- 
cle of 1" radius. No detectable flux is seen at the immediate 
location of the GRB. We place a 2-cr limit in the observed H- 
band magnitude of AB(H) = 26.0 over a 0.5" diameter aper- 
ture for the host galaxy. At z = 3.344, the observed //-band 
magnitude limit allows us to derive a limiting rest-frame ab- 
solute magnitude of M AB (3825)-5 log h > -18.9 for the GRB 
host galaxy. 

4.9. GRB 050922C at z GRB = 2.199 

This burst was detected by Swift (Norris et al. 2005) and 
the optical afterglow was reported by Rykoff et al. (2005) ten 
minutes after the trigger to have R rj 16.0 mag. Spectroscopic 
follow-up was carried out roughly 1 .5 hours after the trigger 
by Jakobsson et al. (2006), who reported z G rb = 2.198 based 
on a series of absorption features. A foreground damped Lya 
absorber (DLA) was also found at z = 2.07 along this sightline. 
Additional high-resolution (8 v w 7 km s" 1 ) echelle spectra of 
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TABLE 2 

A Summary of the Observed Isophotal Magnitudes of Objects at A9 < 1.5" from GRB 050820A 



Objects A6»(") AB(F625W) AB(F775W) AB(F850LP) AB(H) 



Host ... 26.24 ±0.28 26.30 ±0.06 26.04 ±0.05 25.3 ±0.3 

A 1.34±0.05 26.08±0.06 26.11 ±0.05 25.81 ±0.05 > 26.0 

B 0.44±0.05 26.35±0.23 26.21 ±0.05 26.25±0.05 > 26.0 




FIG. 13.— A smoothed H image of the field around GRB050922C at 
Zgrb = 2. 199. No detectable flux is seen at the location of the OT. We place a 
2-cr //-band limiting magnitude of AB(H) = 26 over a 0.5" diameter aperture 
for the host galaxy. 

the afterglow were obtained roughly 3.5 hours after the trig- 
ger using UVES on the VLT Kueyen telescope (D'Elia et al. 
2005). The spectra were retrieved from the ESO data archive 
and processed using our own reduction software. 

The combined echelle spectrum covers a spectral range over 
A = 3300- 10,000 A, allowing accurate estimates of chem- 
ical abudances in the ISM of the GRB host galaxy. Based 
on multiple absorption features from both ground-state and 
excited-state ions, we determine a total neutral hydrogen col- 
umn density of log A^(H I) = 2 1 .5 ± 0. 1 . An absorption-line 
analysis of various ions shows that [S/H] = -2.03 ±0.15, 
[Zn/H] = -2.3 ±0.3 and [Fe/H] = -2.6 ±0.1 (Prochaska et 
al. 2007b; Piranomonte et al. 2008). These measurements to- 
gether imply a dust-to-gas ratio roughly 1 /20 of what is seen 
in the SMC. Adopting the SMC dust-to-gas ratio (Gordon et 
al. 2003), we estimate a visual extinction in the host ISM of 
Ay ~ 0.01 (see also Prochaska et al. 2007b). Searches for H2 
absorption features in the echelle spectrum has also yielded 
null results, placing a 4-cr upper limit on the ISM molecular 
fraction of the host at /h 2 < 10" 6 8 (Tumlinson et al. 2007). 
Additional strong absorbers are found at z = 2.077, z = 2.01 
and z = 1 .99 (Piranomonte et al. 2008). 

We have observed the field around GRB 050922C using 
PANIC on Magellan in August 2007, and obtained a total in- 
tegration of 209 minutes. The mean FWHM of the PSF was 
found to be « 0.5". The final stacked image is presented in 
Figure 13, which has been smoothed using a Gaussian ker- 
nel of FWHM = 0.5". The position of the GRB is marked 
by a circle of 2" radius. No detectable flux is seen at the 
immediate location of the GRB. We place a 2-cr limit in the 
observed //-band magnitude of AB(H) = 26.0 over a 0.5" di- 
ameter aperture for the host galaxy. At z = 2.199, observed 
//-band magnitude limit allows us to derive a limiting rest- 
frame absolute magnitude of Mab{B)-5 log /z > —18.2 for the 
GRB host galaxy. 




FIG. 14.— HST ACS/F814W image of the field around GRB 060206 at 
zgrb = 4.048. The host is detected with AB(F8 14W) = 27.6 ± 0. 1 . Object A 
at AO ss 1" southwest of the host has AB(F814W) = 26.22 ±0.05 and was 
initially noted by Thone et al. (2008a) as a possible host candidate. Given its 
offset from the OT position, we consider galaxy A a foreground galaxy likely 
associated with one of the Mg II absorbers at z = 1 .48 or z = 2.26 found in the 
afterglow spectra. 

Optical images of the field around GRB 050922c have also 
been obtained using LRIS and the g and R c filters on the Keck 
I telescope in July 2006. The mean FWHM of the PSF was 
found to be « 1 .2" in the combined g image and « 1 .0" in the 
combined R c image. No detectable flux is seen at the immedi- 
ate location of the GRB. We place 2-cr limits of AB(g) = 21.0 
and AB(R C ) = 26.4 over a 2" diameter aperture for the host 
galaxy. At z = 2.199, observed g-band magnitude limit al- 
lows us to derive a limiting rest-frame absolute magnitude of 
M Aij (1500)-5 log h > -17.2 for the GRB host galaxy. This is 
comparable to the flux limit in the //-band imaging data based 
on a mean color of UV - B = 1.22 for z = 2 - 3 star-forming 
galaxies (Shapley et al. 2005). 

4. 10. GRB 060206 at zgrb = 4.048 

This burst was detected by Swift (Morris et al. 2006). An 
optical transient was nearly instantaneously identified with 
V w 16.7 (Fynbo et al. 2006b; Boyd et al. 2006). Spectro- 
scopic follow-up of the afterglow was carried out by multiple 
groups (Fynbo et al. 2006a; Prochaska et al. 2006; Aoki et al. 
2006; Hao et al. 2007). The host galaxy at z = 4.048 is found 
to have log N(U I) = 20.85 ±0.10 and [S/H] = -0.85 ±0.10 
based on moderate resolution (5v « 40 kms" 1 ) afterglow 
spectra (Fynbo et al. 2006a; Thone et al. 2008). No H 2 
is detected to a 4-cr limit of /h 2 < 10~ 3,6 (c.f. Fynbo et al. 
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FIG. 15. — A smoothed H image of the field around GRB 060607 at zgrb = 
3.075. No detectable flux is seen at the location of the OT. We place a 2-cr 
//-band limiting magnitude of AB{H) = 26.5 over a 0.5" diameter aperture 
for the host galaxy. 

2006; Tumlinson et al. 2007). In addition to the GRB host, 
two strong Mg II absorbers are found along the line of sight 
at z = 1.48 and z = 2.26 with W(2796) = 0.95 ±0.1 A and 
W(2796) = 1.5 ±0.1 A, respectively (Aoki et al. 2006; Hao et 
al. 2007; Thtine et al. 2008). 

We obtained late-time images of the field around 
GRB 060206 on November 25, 2006, using ACS and the 
F814W filter on board HST (PID=10817). The images were 
processed and registered using the standard pipeline tech- 
nique. A stacked ACS/F814W image is presented in Figure 
14. At the position of the OT (marked by a circle of 0.5" ra- 
dius), we clearly detect a faint source of AZ?(F814W) = 27.6 ± 
0.1. We identify the source as the host galaxy of GRB 060206. 
At z = 4.048, the observed F814W magnitude corresponds 
to a rest-frame absolute magnitude of Mab(1600)-5 log h = 
-17.7 ±0.1 for the GRB host galaxy. 

At AO = 0.96 ±0.02" southwest of the host, we note the 
presence of galaxy A with Afi(F814W) = 26.22 ± 0.05. This 
galaxy was also seen in a deep, ground-based /-'-band im- 
age published in Thone et al. (2008). The authors initially 
identified galaxy A as a candidate for the GRB host galaxy, 
but revised the identification after analyzing the HST images. 
Given the large angular distance to the OT (corresponding to 
p « 4.8 h kpc at z = 4.408) and the presence of at least two 
strong Mg II absorbers in the afterglow spectra, we conclude 
that galaxy A is likely a foreground galaxy associated with 
one of the Mg II absorbers. 

4.11. GRB 060607 at zgrb = 3.075 

This burst was detected by Swift and the optical afterglow 
was detected one minute after the trigger by UVOT and the 
white filter (1600-6500 A) on board the satellite (Ziaeepour 
et al. 2006). The initial brightness was estimated ea 15.7 mag 
(Ziaeepour et al. 2006). A series of high-resolution (<5v w 7 
km s" 1 ) optical spectra of the afterglow were obtained using 
UVES on the VLT telescope by Ledoux et al. (2006), start- 
ing 7.5 minutes after the trigger. These authors identified 
the GRB host at z = 3.082 and noted two probable DLAs at 
z = 2.937 and z = 3.05. We retrieved the echelle spectra from 
the ESO data archive and processed the data using our own 
reduction software. 

The afterglow spectrum of GRB 060607 exhibits a number 
of unusual features. First, the GRB host appears to arise in 



a weak Lya absorber at z = 3.075, with associated Si IV and 
C IV absorption feature but no trace of low-ionization species. 
We determine log A^(H I) = 16.85 ±0.10 based on a simulta- 
neous Voigt profile analysis of the Lyman absorption series, 
using the VPFIT 3 software package. This system represents 
the only optically thin absorber detected in a GRB host galaxy 
(Jakobsson et al. 2006; Chen et al. 2007a). Second, the strong 
Lya absorber at z = 3.05, « 1840 km s" 1 blueshifted from the 
host, is a Lyman limit system of logAf(H I) = 19.2 ±0.1 with 
weak metal absorption lines detected, implying a metallicity 
< 1/100 solar. Third, the strong Lya absorber at z = 2.937 ex- 
hibits complex kinematic profiles in various ionic transitions, 
spreading over Av ss 600 km s _1 . We measure log N(H I) = 
19.6±0.1. 

We have observed the field around GRB 060607 using 
PANIC on Magellan in August 2007, and obtained a total in- 
tegration of 309 minutes. The mean FWHM of the PSF was 
found to be w 0.5". The final stacked image is presented in 
Figure 15, which has been smoothed using a Gaussian kernel 
of FWHM = 0.5". The position of the GRB is marked by a 
circle of 2" radius. Despite the discovery of three strong Lya 
absorbers at z = 2.93-3.08, no detectable flux is seen at the 
immediate location of the GRB. We place a 2-cr limit in the 
observed H-band magnitude of AB(H) = 26.5 over a 0.5" di- 
ameter aperture for the host galaxy. At z = 3.075, observed 
H-band magnitude limit allows us to derive a limiting rest- 
frame absolute magnitude of Mab{B)-5 log > —18.3 for the 
GRB host galaxy. 

Optical images of the field around GRB 060607 have also 
been obtained using LRIS and the g and R c filters on the Keck 
I telescope in July 2006. The mean FWHM of the PSF was 
found to be w 1.0" in the combined g image and w 1.5" 
in the combined R c image. The effective seeing was com- 
promised in the R c images due to a focus problem with the 
LRIS red-side. No detectable flux is found at the immedi- 
ate location of the GRB. We place 2-cr limits of AB(g) = 26.8 
and AB(R C ) = 26.5 over a 2" diameter aperture for the host 
galaxy. At z = 3.075, observed /? c -band magnitude limit al- 
lows us to derive a limiting rest-frame absolute magnitude of 
M AJJ (1600)-5 log h > -18.3 for the GRB host galaxy. 

4.12. GRB 07072 IB at z G rb = 3.626 

This burst was detected by Swift (Ziaeepour et al. 2007) and 
the optical afterglow was detected one minute after the trig- 
ger by UVOT and the white filter (1600-6500 A) on board 
the satellite (Schady et al. 2007). The initial brightness was 
estimated m 15.9 mag (Schady et al. 2007). Spectroscopic 
follow-up was carried out by Malesani et al. (2007), who re- 
ported zgrb = 3.626 based on the presence of a strong Lya 
feature and a series of metal absorption lines. The host galaxy 
is found to have logA^(HI) = 21.50 ±0.20 based on mod- 
erate resolution afterglow spectra (Jakobsson et al. 2008, in 
preparation). A additional intervening damped Lya absorber 
is found at z = 3.09. 

We have observed the field around GRB 07072 IB using 
PANIC on Magellan in August 2007, and obtained a total in- 
tegration of 172 minutes. The mean FWHM of the PSF was 
found to be w 0.5". The final stacked image is presented in 
Figure 16, which has been smoothed using a Gaussian ker- 
nel of FWHM = 0.5". The position of the GRB is marked by 
the cross within a circle of 2" radius. No detectable flux is 
seen at the immediate location of the GRB. However, an ex- 

3 http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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FIG. 16.— A smoothed # image of the field around GRB070721B at 
ZGRB = 3.626. No detectable flux is seen at the location of the OT (marked by 
the cross). An extended source is seen at AS = 0.9" sourtheast of the OT with 
AB(H) = 23.7 ±0.1, which has been confirmed to be at z = 3.09 and associ- 
ated with the foreground DLA along the line of sight (Fynbo 2008, private 
communication; Milvang- Jensen et al. 2008, in preparation). We place a 2-<r 
ff-band limiting magnitude of AB(H) = 25.8 over a 0.5" diameter aperture 
for the host galaxy. 




FIG. 17.— Optical images of the field around GRB000301C at z G rb = 
2.0404. The epoch during which the images were taken is indicated at the 
bottom of each panel. The late-epoch image (right panel) has been smoothed 
using a Gaussian kernel of FWHM = 0.15", which is roughly the size of the 
PSF. At the location of the OT, we detect faint emission of the host galaxy in 
the late-time image and estimate AB(clear) = 28.9 ± 0.5 over a 0.5" diameter 
aperture (consistent with what was reported in Fruchter et al. 2006). 

tended source is seen at A9 = 0.9" sourtheast of the OT with 
AB(H) = 23.7 ±0.1. This source has been confirmed to be at 
z = 3 .09, associated with the foreground DLA along the line of 
sight (Fynbo 2008, private communication; Milvang-Jensen 
et al. 2008, in preparation). This represents by far the most lu- 
minous DLA galaxy found at z > 2 (e.g. M0ller et al. 2002a). 
At z = 3.09, the projected distance between the foreground 
DLA and the extended source would be p = 4.8 h kpc. We 
therefore consider the GRB host galaxy missing in our H- 
band image and place a 2-er limit of AB(H) = 25.8 over a 0.5" 
diameter aperture for the host galaxy. At z = 3.626, observed 
//-band magnitude limit allows us to derive a limiting rest- 
frame absolute magnitude of M,4b(3600)-5 log h > -19.3 for 
the GRB host galaxy. 

4.13. GRBs 000301C, 000926, & 021004 

Early-time imaging and spectroscopic observations of the 
afterglow of GRB 00030 1C have been presented in Jensen et 
al. (2001), who reported based on low-resolution afterglow 
spectra that the source is at z = 2.0404. The host ISM in front 



of the afterglow was found to have logA^(H I) = 21.2±0.5. 
Multi-epoch imaging follow-up was carried out using STIS 
on board HST and the clear filter under program #8189 (PI: 
Fruchter; see also Fruchter et al. 2006). We retrieved available 
imaging data from the HST data archive and analyzed the im- 
ages ourselves. In the left panel of Figure 17, we present a 
stacked image of the field obtained in April 2000 with a to- 
tal exposure time of 9391 s. The optical transient is clearly 
visible in this early epoch image. In contrast, the images (to- 
tal exposure time of 7031 s) obtained in February 2001 (right 
panel of Figure 17) exhibit only faint emission at the location 
of the GRB. Attributing the faint emission to the host galaxy, 
we measure a total flux of AB(clear) = 28.9 ± 0.5 over a 0.5" 
diameter aperture. This is consistent with the measurement 
reported by Fruchter et al. (2006). Accounting for the band- 
pass difference (e.g. Figure 3), we derive AB{R) = 28.8 ±0.5 
and M AB (2200)-5 log h = -15.2 ±0.5. 

Both GRB 000926 and GRB 021004 have been studied ex- 
tensively by previous authors. Their host galaxies are identi- 
fied in early effort to carry out imaging follow-up with HST. 
Here we briefly review known emission and absorption prop- 
erties of the two host galaxies to complete the discussion of 
individual GRB hosts. 

Detailed spectroscopic and imaging studies of GRB 000926 
are presented in Harrison et al. (2001), Fynbo et al. (2001), 
and Castro et al. (2003). The afterglow spectrum displays 
a DLA of log/Y(H I) = 21.3 ±0.3 and a series of metal ab- 
sorption features at z = 2.0385 with an estimated metallic - 
ity of [Zn/H] =-0.17 ±0.15 and metal abundance ratio of 
[Zn/Fe] = +1 .3 ± 0. 15 (Castro et al. 2003). Adopting the dust- 
to-gas ratio and the SMC extinction law, we derive a visual 
extinction of A v = 0.15 and E(B-V) = 0.06. 

The host galaxy is identified both in Lya emission on 
the ground (Fynbo et al. 2002) and in space images ob- 
tained using HST WFPC2 and the F606W filter (Castro et 
al. 2003). It displays a disturbed morphology at rest-frame 
UV wavelengths with extended emission over a w 15 h~ l 
kpc projected size. Additional near-infrared /-band images 
show that the host has AB(J) = 24.1^ 4 (Christensen et al. 
2004) which, together with available WFPC2 photometry, 
constrains its absolute magnitudes at rest-frame UV and opti- 
cal wavelengths. Including the intrinsic extinction correction 
estimated from absorption line analysis, we find Mab(2000)- 
5 log h = -19.63 ± 0.07 and M AB (B) - 5 log h = -20.26^. 
Comparing the observed SED with a suite of stellar synthetic 
models described in § 4.6, we find the host galaxy is best de- 
cribed with a mean stellar age of ss 570 Myr and an exponen- 
tially declined SFR of e-folding time 1 Gyr (left panels of Fig- 
ure 18). Adopting the best-fit stellar synthetic model, we can 
constrain the total underlying stellar mass to be = 2.1^" x 
10 9 h~ 2 M Q , consistent with M* = (1.6 ± 3.1) x 10 9 h~ 2 M Q 
estimated by Savaglio et al. (2008). 

Detailed spectroscopic and imaging studies of GRB 021004 
are presented in M0ller et al. 2002b, Mirabal et al. 2003, 
Schaefer et al. 2003, Fiore et al. 2005, and Fynbo et al. 2005. 
The afterglow spectrum displays a combination of emission 
and absorption due to the hydrogen Lya transition and a series 
of metal absorption lines at z = 2.329. The presence of Lya 
emission in the afterglow spectrum makes it difficult to deter- 
mine A^(H I) precisely. Based on the absence of Lyman limit 
photons, Fynbo et al. (2005) derive log N(H I) = 19.5 ± 0.5. 
This is therefore one of the four GRB Lyman limit absorbers 
published so far (including GRB 060526 in Jakobsson et al. 
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2006). We have retrieved available echelle spectra of the OT 
from the ESO data archive. Based on the absorption strengths 
of Si II and Fe II observed in the combined spectrum, we take 
into account necessary corrections for the ionization fraction 
of the gas and place an upper limit to the metallicity of the gas 
at[Si/H]<-l. 

Optical and near-infrared photometry of the host galaxy 
of GRB 021004 are reported and analyzed in Fynbo et al. 
(2005). The authors estimate a mean stellar age of « 42 
Myr. Adopting the optical and near-infrared photometric 
measurements of Fynbo et al. (2005), we derive a rest-frame 
absolute magnitude of M AB (B) - 5 log h = -20.38 ±0.15 or 
M AB (2000)-5 log h = -19.83 ±0.07. Comparing the observed 
SED with a suite of stellar synthetic models described in § 4.6, 
we find the host galaxy is best decribed by a single starburst 
episode that occurred w 39 Myr ago, and constrain the total 
underlying stellar mass to be M* = 1 .3^ 5 x 10 9 h~ 2 M Q . This 
is smaller than M* = (7.8 ±0.7) x 10 9 h~ 2 M Q estimated by 
Savaglio et al. (2008). The discrepancy may be understood 
by the poorly constrained star formation history, displayed in 
the right panels of Figure 18. 

We note that two strong Mg II absorbers are found at z = 
1.38 and z = 1.60 with W(2796) = 1.6 A and W(2796) = 1.4 
A, respectively. The large absorption width of these fore- 
ground absorbers imply possible contaminating light in the 
host imaging observation (c.f. Pollack et al. 2008). In- 
deed, the host galaxy is clearly resolved into two compact 
clumps in available HST images (e.g. Fynbo et al. 2005; 
Chen et al. 2007b), but extended Lya emission is detected 
in ground-based narrow-band imaging follow-up (Jakobsson 
et al. 2005). In the following discussion, we consider the de- 
rived absolute magnitude as an upper limit to the intrinsic lu- 
minosity of this host galaxy. 

5. DISCUSSION 

We have carried out an optical and near-infrared imaging 
survey of the fields around 15 GRBs at z > 2. The GRBs 
are selected to have early-time afterglow spectra in order to 
compare ISM absorption-line properties with stellar proper- 
ties. The redshifts of the GRBs span a range from z = 2.04 
to z = 4.05, and the neutral hydrogen column densities of 
the GRB host ISM span a range from log/V(H I) = 16.9 to 
log N(U I) = 22.6 (Figure 1). 

In addition to the five previously known GRB host galaxies, 
we report new detections for the host galaxies of GRB 050820 
and GRB 060206. The seven identified GRB host galax- 
ies have rest-frame UV absolute magnitudes spanning from 
M AB (U)-5 log h = -15.2 to M AB (UV)-5 log h = -19.8 mag. 
For sources with a known dust-to-gas ratio from afterglow 
absorption-line analysis, we correct the rest-frame UV mag- 
nitude for dust extinction assuming an SMC extinction law. 
The inferred SFR spans a range from 0.6 to 3.8 h~ 2 M 
yr" 1 . We are able to constrain the underlying stellar popula- 
tions for three host galaxies, GRB 000926, GRB 021004, and 
GRB 050820, based on comparisons of the observed optical 
and near-infrared broad-band colors and a suite of stellar pop- 
ulation synthetic models . We estimate total stellar masses of 
between M» = (0.9-2.1) x 10 9 h~ 2 M Q for these hosts. We 
also place 2-a upper limits for the rest-frame luminosities of 
the remaining eight GRB host galaxies based on the depths 
in available optical and near-infrared images. Finally, high 
spatial resolution images from the HST allow us to deblend 
the GRB host galaxies from foreground absorbers and to un- 
veil a range of rest-frame UV morphology between compact 
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FIG. 18. — Constraints on the galaxy age and star formation history based 
on a stellar synthesis model analysis. Left: The broad-band spectral energy 
distribution of the host galaxy of GRB 000926 at z = 2.0379 (points with 
eiTorbars), in comparison to the best-fit template (solid histogram) and the 
predicted broad-band photometric points (open squares). The best-fit SED is 
charaterized by a declining SFR of e-folding time r = 300 Myr, solar metal- 
licity, and intrinsic dust extinction of E(B- V) = 0.025 following the SMC 
extinction law. The likelihood function of the stellar age for the host galaxy 
is display in the bottom left panel, suggesting that a recent starburst occurred 
about 650 Myr prior to the GRB explosion. Right: A summary of the stellar 
synthesis model analysis for the host galaxy of GRB 021004 at z = 2.329. 
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(e.g. GRB 060206) and extended (e.g. GRB 050820) emission 
features of the hosts. 

A summary of known absorption-line properties and stellar 
properties of the 15 GRB host galaxies is presented in Table 
3. Combining early-time, high-resolution afterglow spectra 
and the results of late-time imaging survey of the GRB fields 
allows us to address a number of issues regarding both the na- 
ture of GRB progenitor environment and star-forming physics 
in distant starburst galaxies. 

5.1. The Luminosity Distribution of GRB Host Galaxies 

First, we examine the luminosity distribution of GRB host 
galaxies above z = 2 based on the survey result of our sam- 
ple. The goal is to characterize the nature of GRB host galax- 
ies based on comparisons of their luminosity distribution and 
the luminosity distribution of a randomly selected sample 
from the field galaxy population. We incorporate the non- 
detections in our follow-up imaging survey by evaluating the 
cumulative maximum fraction J- mdx of GRB host galaxies that 
are fainter than a given UV absolute magnitude M max (t/V) 
(solid histogram in Figure 19). For the host of GRB 030429, 
we have a constraint only for the rest-frame Z?-band magni- 
tude. We infer its corresponding UV magnitude based on the 
mean color of (UV-B) = 1.22 (with an r.m.s. scatter of 0.3 
mag) observed for luminous starburst galaxies at z = 2 - 3 in 
Shapley et al. (2005). For the host of GRB 060607, we have 
estimated limiting magnitudes both in the rest-frame UV and 
B bands. We adopt the more sensitive limit based on the con- 
version of UV-/3 = 1.22 4 . The empirical observations of the 
sample of 15 GRB fields confirm the previous understanding 
for z ~ 1 GRB hosts (e.g. Le Floc'h et al. 2003): as much as 
70% of long-duration GRBs may originate in galaxies fainter 
than 0.1 L*. 

At z = 2 - 4, the field galaxy population is now well char- 
acterized by a Schechter luminosity function (f>(M) at rest- 
frame UV wavelengths with Mab* - 5 log h = -20.2 ± 0.1, 
(f>* = (4 ±0.6) x 10" 3 /r 3 Mpc~ 3 , and a faint-end slope a = 
-1.7 ±0.1 (Bouwens et al. 2007; Reddy et al. 2008). If the 
GRB host galaxies are representative of the field galaxy pop- 
ulation, then we expect that the fraction of host galaxies found 
in a luminosity interval is proportional to the space density of 
galaxies in the luminosity range. The expected cumulative 
maximum fraction of the host galaxies versus UV magnitude 
can be estimated following 



^max[M>M max ] OC 
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where the first term extends over n known host galaxies 
that have Mab(UV) > M max , the second term extends over m 
unidentified host galaxies that may be as luminous as M max , 
and 4>{M) is the galaxy luminosity function. The luminos- 
ity distribution is expected to resemble the galaxy luminos- 
ity function with a dominant fraction attributed to faint dwarf 
galaxies (e.g. Jakobsson et al. 2005; Fynbo et al. 2008). 

We test this hypothesis by calculating J- m - dx for a sample 
of 15 random galaxies that share the known luminosities of 

4 We note that while there is no apparent trend between UV —B and Muv 
in the luminous starburst sample of Shapley et al. (2005), it is possible that 
that fainter dwarf starburst galaxies may be bluer. Given that only two fields 
(GRB 030429 and GRB 060607) are affected by this conversion, we find that 
the results presented in this section are not sensitive to the adopted UV — B 
color. 
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FIG. 19. — Cumulative distribution of the maximum fraction of GRB 
host galaxies that are fainter than a given UV absolute magnitude (solid his- 
togram), based on the seven identified host galaxies and eight fields with up- 
per limits in our imaging sample. The curves show model expectations based 
on different hypothesis for the origin of GRB host galaxies. We consider 
two scenarios: (1) GRB host galaxies are representative of the general galaxy 
population and (2) GRB host galaxies originate preferencially in galaxies of 
higher star formation rate. Adopting the best-fit UV luminosity function of 
star-forming galaxies at z = 3 - 4 from Bouwens et al. (2007) and Reddy 
et al. (2008), we show the expectations for scenario (1) in dashed curves 
for two different minimum luminosity cutoffs, L mm = 0.01 L* (bottom) and 
Lmin = 0.005 L» (top). The expectation for scenario (2) with L mm = 0.001 L* 
is shown in the solid curve, and the dotted curves are for L m ; n = 0.01 L* 
(lower) and L m j n = 0.0001 L» (upper). The comparisons provide empirical 
evidence supporting the expectation that GRB host galaxies form an SFR 
weighted field galaxy sample. 

the seven identified GRB host galaxies and the empirical 
2-cr upper limits for the remaining eight GRB host galax- 
ies. We experiment with different M m \ n (with corresponding 
minimum luminosity L mm ) in Equation (1). The results for 
Lmin = . 1 L* and for L m j n = . 005 L* are shown as the dashed 
curves in Figure 19. Smaller L m \ n models predict a larger con- 
tributions from fainter galaxies. The model expections for dif- 
ferent adopted L m ; n clearly deviate from the observed distri- 
bution. The hypothesis that the host galaxies of long-duration 
GRBs trace random field galaxies is rejected at > 99% confi- 
dence level based on a Kolmogorov-Smirnov (KS) test. 

Next, under the assumption that GRBs trace instantaneous 
star formation, galaxies with higher on-going SFR are ex- 
pected to have higher probability to host a GRB event. Ap- 
plying the rest-frame UV luminosity as a measure of the on- 
going SFR, we modify Equation (1) to include a UV lumi- 
nosity weighting, J L(f)(L)d(L), for calculating the expected 
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The results are shown as solid and dotted curves in Figure 19. 
The observations are best described under this hypothesis for 
L m in = 0.001 L* (solid curve) based on a KS test, but we cannot 
rule out models with L m ; n = 0.01 L* (lower dotted curve) or 
Lmm = 0.0001 L* (upper dotted curve). 

In summary, the sample of seven known GRB host galax- 
ies and eight upper limits for unidentified hosts allows us to 
determine the cumulative maximum fraction of the z > 2 host 
galaxy population as a function of rest-frame UV magnitude. 
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Adopting rest-frame UV luminosity as a measure of on-going 
SFR, we find that the empirical sample is best described by a 
SFR- weighted sample of the field galaxy population. Models 
that do not include SFR weighting can be ruled out at > 99% 
confidence level. Based on the best-fit SFR-weighted model, 
we estimate a median luminosity for the GRB host galaxies at 
«0.1L*. 

A similar analysis has been presented in Jakobsson et al. 
(2005), who derived constraints for the luminosity function 
of GRB host galaxies under the assumption that the ob- 
served brightness distribution of the hosts follows a luminos- 
ity weighted field galaxy population. Our study differs from 
the approach of Jakobsson et al. in that our analysis is based 
on a uniform set of photometric measurements from our own 
imaging survey. Then we adopt the known UV luminosity 
function of the field galaxy population (with a faint-end slope 
of a = -1 .7) and examine different hypotheses for generating 
the observed GRB host galaxy sample. We confirm that the 
host galaxy population is representative of a UV luminosity 
weighted sample. The difference between the host luminosity 
function of Jakobsson et al. and the best-fit luminosity of star- 
forming galaxies at z > 2 may be due to uncertainties in pub- 
lished photometric data of the host galaxies in their sample 
and uncertainties in the line-of-sight absorption properties. 

We note that although the selection criterion of our GRB 
sample is based on available afterglow spectra that presum- 
ably includes only GRBs with relatively bright optical after- 
glows, the broad range in the isotropic energy release of the 
GRBs (see Table 1 in Chen et al. 2007a) indicates that the 
GRBs in our sample are not an overly biased portion of the 
long-duration GRB population. That is, we have selected a 
representative subsample of the unobscured GRB host galaxy 
population. However, the presence of 20% dark bursts that do 
not have optical afterglows found (e.g. Tanvir & Jakobsson 
2007) suggests that some GRBs originate in heavily obscured 
star-forming regions that are not included in our sample and 
are likely missed in the UV selected field galaxy sample as 
well. To constrain the fraction of dust obscured star-forming 
galaxies requires an independent study of the fields around 
dark bursts. 

5.2. The Luminosity-Metallicity Relation in GRB Host 
Galaxies 

We combine known absorption-line properties with esti- 
mated galaxy UV luminosity to investigate the gas kinematics 
and chemical enrichment in the ISM of GRB host galaxies. 
The goals of this study are (1) to examine the physical pro- 
cesses that determine the observed absorption-properties in 
the host ISM, (2) to investigate whether there is a metallicity 
cutoff in GRB host galaxies as favored by various theoreti- 
cal models (e.g. Hirschi et al. 2005; Yoon & Langer 2005; 
Woosley & Heger 2006), and (3) to probe the luminosity- 
metallicity relation below the magnitude limit of most previ- 
ous studies at z > 2 (e.g. Erb et al. 2006; Maiolino et al. 2008). 

We first compare rest-frame absorption equivalent widths of 
the Sill A 1526 transition W(1526) observed in the host with 
the absolute UV magnitude. For this study, we consider only 
host DLAs because the lyman limit absorbers do not include 
neutral ISM in the host galaxies. The Si II A 1526 transition is 
typically saturated at W(1526) > 0.3 A. The line width pro- 
vides a measure of the velocity field of cool clouds along the 
line of sight through a galactic halo, rather than the total gas 
column density (e.g. Prochaska et al. 2008). We adopt the 
W(1526) measurements published in Prochaska et al. (2008). 
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FIG. 20. — Absolute UV magnitude of GRB host galaxies versus rest-frame 
absorption equivalent width of the Si II A 1526 transition in the hosts as ob- 
served in afterglow spectra. For this study, we consider only host DLAs be- 
cause the lyman limit absorbers do not include neutral ISM in the host galax- 
ies. The Sill A 1526 transition is typically saturated at rest-frame absorption 
equivalent width IV(1526) > 0.3 A. The line width therefore provides a mea- 
sure of the velocity field of cool clouds along the line of sight through a halo 
(e.g. Prochaska et al. 2008). We find a moderate trend with stronger W(1526) 
appearing in more luminous hosts at nearly 95% confidence level. 

Including additional measurement for GRB 060206 in public 
afterglow spectra from Subaru Science Data Archive (Aoki et 
al. 2006), we have assembled eight GRB hosts galaxies for 
this study. 

Figure 20 shows that there exists a relatively significant 
trend with stronger Si II transitions appearing in more lumi- 
nous hosts. Including upper limits, we find based on a gener- 
alized Kendall test that the probability of a positive correlation 
between M uv and W(1526) is nearly 95%. While the posi- 
tive correlation is similar to the mass-metallicity relation seen 
in field galaxies (e.g. Erb et al. 2006; Maiolino et al. 2008) 
, we note that half of the host galaxies have W(1526) > 1.5 
A, corresponding to a velocity width of ~ 300 km s" 1 . At- 
tributing the large line width to gravitational motion of clouds 
within the host dark matter halos would require a halo mass 
> 10 11 Ii~ 1 Mq, comparable to the mass scale found for ha- 
los that host luminous starburst galaxies at z ~ 2 (e.g. Conroy 
et al. 2008) but at odds with expectations for starburst galax- 
ies of low stellar mass and low luminosity. But because M uv 
serves as a measure of the on-going SFR, the observed M uv 
versus W(1526) correlation implies that the velocity field ob- 
served in GRB host galaxies is driven by galactic outflows. 

Next, we examine whether there exists a correlation be- 
tween host luminosity and ISM metallicity, in comparison 
to what is known for luminous starburst galaxies published 
by Erb et al. (2006). We note two apprarent caveats in our 
study. First, chemical abundances of the GRB host galax- 
ies in our sample are determined for the cold neutral medium 
using absorption-line techniques, whereas the metallicities of 
field galaxies are determined from the integrated emission- 
line fluxes of their HII regions. Although extensive stud- 
ies have yet to be undertaken to examine possible systematic 
differences between absorption- and emission-line abundance 
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measurements, available evidence based on limited studies of 
nearby starburst galaxies have yielded consistent metallicity 
measurements using either absorption-line or emission-line 
techniques (see Russell & Dopita 1992 and Welty et al. 1997, 
1999 for the Large and Small Magellanic Clouds; Lecavelier 
des Estangs et al. 2004 for IZw 18; and Schulte-Ladbeck et 
al. 2005 and Bowen et al. 2005 for SBS 1543+593). In the 
few cases where emission and absorption line measurements 
have been compared in galaxies at z > 2, the two methods 
have been found to give consistent answers to within a factor 
of ~ 2 (e.g. Pettini 2006). 

Second, absorption-line measurements are for the ISM 
along the afterglow line of sight, whereas emission-line ob- 
servations are integrated measurements averaged over the en- 
tire galaxies. A large metallicity gradient is commonly seen in 
nearby galaxies with a slope varying between dZ/dr w -0.02 
dex / kpc to -0.07 dex / kpc (e.g. Zaritsky et al. 1994; van Zee 
et al. 1998; Kennicutt et al. 2003). Therefore, a metallicity 
measurement for individual sightlines may not be representa- 
tive of the global mean value of the entire host galaxy. We 
note, however, that z > 2 galaxies are relatively compact with 
typical half light radii of w 2 kpc (Bouwens et al. 2004; Law 
et al. 2007). Previous observations (e.g. Bloom et al. 2002; 
Fruchter et al. 2006) and current findings for GRB 050820 and 
GRB 060206 show that the GRBs occur within < 2 kpc radius 
of their host galaxies. This is also consistent with the large 
N(H I) observed in the afterglow spectra with the exception of 
four non-DLAs. Spatial variation of the observed metallicity 
across distant star-forming galaxies is therefore not expected 
to exceed 0.2 dex. We proceed with an analysis that compares 
the ISM metal content of GRB galaxies with those of known 
starburst galaxies at z ~ 2. 

Figure 21 shows the the luminosity-metallicity relation re- 
produced from Erb et al. (2006) for starburst galaxies at 
(z) = 2.3 (open squares). The oxygen abundances are mea- 
sured using the N2 index, which is a measure of the flux ratio 
between the observed [Nil] A6584 to Ha lines 5 . The rest- 
frame B-band magnitudes published in Erb et al. have been 
converted to UV magnitude according to UV = B + 1 .22 to 
facilitate comparisons with the GRB host galaxy population 
and with predictions from numerical simulations (crosses in 
Figure 21, reproduced from Fynbo et al. 2008). The UV -B 
color conversion is estimated based on the mean R-K colors 
for these galaxies in Shapley et al. (2005). 

We include in Figure 21 measurements for nine GRB host 
galaxies for which measurments (solid points) or constraints 
(open circles) on the ISM metallicity are available. Only 
known host DLAs are considered in the comparison here. Left 
arrows represent the fields (GRB 050401, GRB 050922c and 
GRB 050730) for which no hosts have been found and 2-cr 
upper limits to the intrinsic luminosity are shown. The lower 
limits represent those GRB fields for which only moderate 
resolution afterglow spectra are available and the observed 
absorption-line strengths represent a lower limit to the intrin- 
sic values (e.g. Prochaska 2006). 

Three interesting features are seen in Figure 21. First, the 
majority of the GRB host galaxies are fainter than the faintest 
starburst galaxies studied in magnitude-limited surveys (open 
squares) and offer a promising probe to extend the studies of 

5 See Erb et al. 2006 for a discussion of possible systematic uncertainties 
associated with the N2 index. While it is known that the N2 index saturates at 
solar metallicity, the saturation does not affect our comparison because GRB 
host galaxies appear to have mostly sub-solar abundances (e.g. Prochaska et 
al. 2007a). 
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FIG. 21. — The luminosity-metallicity relation in distant galaxies. Solid 
points with errorbars are GRB host galaxies with accurate absorption-line 
metallicity measurements, while open circles are GRB host galaxies with 
constraint on their ISM metallicities due to either insufficient spectral reso- 
lution in the afterglow spectra or unknown ionization fraction of the ISM. 
Open squares are emission-line metallicity measurements for luminous star- 
burst galaxies at z ~ 2 from Erb et al. (2006). The oxygen abundance was 
evaluated using the observed [N II] to Ha line ratio. We have converted the 
rest-frame B-band magnitude to UV magnitude according to UV = B + 1.22, 
which is estimated based on the mean R—K colors for these galaxies in Shap- 
ley et al. (2005). Crosses represent the predicted luminosity-metallicity re- 
lation for starburst galaxies at z. ~ 3 (Fynbo et al. 2008) based on numerical 
simulations by Sommer-Larsen & Fynbo (2008). 

ISM metal enrichment to fainter luminosity limits than exist- 
ing faint galaxy surveys (Djorgovski et al. 2004). Second, the 
metallicities found in GRB host galaxies span two orders of 
magnitude from -2.2 dex below solar to ks -0.2 below solar, 
while the rest-frame UV luminosity of the galaxies extends 
from sa 1/2L* to fainter than 0.05 L». There is no apparent 
upper metallicity cutoff in the sample of GRB host galax- 
ies at z > 2. Finally, we include predictions from numerical 
simulations performed by Sommer-Larsen & Fynbo (2008) 
that include SNe feedback. The crosses are reproduced from 
Fynbo et al. (2008) for comparison with observations. De- 
spite the presence of lower limits in the metallicities of three 
GRB hosts, there is a moderate trend that suggests a steeper 
luminosity-metallicity relation than what is seen in simula- 
tions. We note that this is unlikely due to a selection bias 
against metal enriched faint galaxies based on the study pre- 
sented in § 5.1, which shows that the GRB host galaxies are 
consistent with an SFR selected field galaxy sample. 

Different theoretical studies have been carried out to un- 
derstand the key astrophysical processes that determine the 
luminosity-metallicity relation found in galaxies from differ- 
ent epochs (e.g. Dekel & Woo 2003; Tassis et al. 2008; Brooks 
et al. 2007; Finlator & Dave 2008). The steep decline in ISM 
metallicity toward fainter magnitudes may be interpreted as a 
signature of SNe driven outflows that removes metals more ef- 
fectively in lower-mass halos than in more massive ones (e.g. 
Dekel & Woo 2003), or as a signature of inefficient star for- 
mation in low mass halos due to low gas density (e.g. Tas- 
sis et al. 2008; Robertson & Kravtsov 2008). The model ex- 
pectations presented in Figure 20 includes SNe feedback but 
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no ISM radiation field to account for subsequent destruction 
of molecules in star-forming regions. The observed steeper 
luminosity-metallicity relation in faint galaxies and the ab- 
sence of molecular gas in GRB host ISM (Tumlinson et al. 
2007) suggests that an enhanced ISM radiation field from 
young stars may be important to effectively reduce subsequent 
star formation and ISM chemical enrichment in dwarf galax- 
ies. 

5.3. Empirical Constraints for the ISM Radiation Field in 

GRB Host Galaxies 

To understand the lack of molecular gas in GRB host ISM, 
Tumlinson et al. (2007) performed a comparison between ob- 
servations and predictions from a grid of models that cover a 
parameter space spanned by four unknown properties. These 
include gas density, metallicity, clouds size and interstellar ra- 
diation field. The authors showed that the lack of molecular 
gas in the host ISM can be understood by a combination of 
low metallicity (low grain production rate) and high interstel- 
lar UV radiation field in the host ISM. Recently, Whalen et al. 
(2008) carried out a set of numerical simulations to examine 
whether the absence of molecular gas is due to the afterglow 
radiation field or an enhanced interstellar UV radiation field 
from a pre-existing H II region where the progenitor star re- 
sides. Taking into account known metallicities from afterglow 
absorption analysis, Whalen et al. concluded that, similar to 
what is found in the Magellanic Clouds, ISM radiation fields 
of intensities up to 100 times the Galactic mean are necessary 
to explain the absence of H2 in the host ISM. 

With the resolved galaxy morphologies seen in the HST 
images of GRB 000926, GRB 030323, GRB 050820A, and 
GRB 060206, we can measure directly the mean radiation 
field in the ISM of the host galaxies at near UV wavelengths 
(« 1500-2000 A). The mean UV radiation intensity is de- 
termined by averaging the observed UV flux over the extent 
of each host galaxy. For comparison with models, we con- 
vert from the observed mean radiation intensity at near UV 
wavelengths to the far UV Lyman- Warner band (11.8-13.6 
eV) based on the Galactic UV radiation spectrum of Gond- 
halekar et al. (1980). The estimated far UV radiation fields Iq 
of individual host galaxies are presented in Table 3. Adopt- 
ing a mean Galactic radiation field at far UV wavelengths of 
/o = 2.5 x 10" 8 photons cm" 2 s" 1 Hz" 1 (Gondhalekar et al. 
1980), we find that, despite a moderate SFR, the ISM ra- 
diation field in the GRB host galaxies spans a range over 
35 - 350 x the Galactic mean value due to a relatively com- 
pact size. The observations provide empirical support for the 
previous theoretical understanding (e.g. Whalen et al. 2008) 
that strong UV radiation fields in low-metallicity ISM envi- 
ronment may be a dominant factor for suppressing the forma- 
tion of molecules and therefore subsequent star formation in 
dwarf galaxies. 

5.4. Properties of Foreground Mg II Absorbing Galaxies 

A puzzling observation of the absorption properties along 
GRB lines of sight is the apparent overdensity of strong (rest- 
frame absorption equivalent width W(2796) > 1 A) Mg II ab- 
sorbers (Prochter et al. 2006). Afterglow spectra exhibit on 
average w 4 x more strong Mg II absorbers than random QSO 
spectra (Prochter et al. 2006; Sudilovsky et al. 2007), al- 
though such over abundance is not seen in CIV absorbers 
at somewhat higher redshift (Sudilovsky et al. 2007; Tejos 
et al. 2007). Various scenarios have been considered to ex- 
plain this discrepancy, including dust extinction due to the 
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FIG. 22. — Cumulative number of galaxies versus angular distance to the 
afterglow line of sight for GRB fields with (top) and without (bottom) inter- 
vening strong Mg II absorbers found in afterglow spectra. All seven fields 
have high spatial resolution HST images available for identifying galaxies 
brighter than AB(R) = 28. In all four fields with known intervening strong 
Mg II galaxies, we find additional objects at A 8 < 2" from the GRB after- 
glow position. In contrast, none is seen at this small angular separation in 
fields without known Mg II absorbers. 

presence of these absorbers that biases observations of QSO 
sightlines, the GRBs being gravitationally lensed by the fore- 
ground absorbers, and the absorbers being intrinsic to the 
GRBs (Prochter et al. 2006; Porciani et al. 2007). However, 
none of these scenarios alone is found sufficient to explain 
the observed overabundance along afterglow sightlines. A 
different scenario has been proposed by Frank et al. (2007), 
who consider different beam sizes between QSOs and GRB 
afterglows as a possible explanation, but two important con- 
sequences are associated with this scenario. First, a partial 
covering of Mg II gas is expected along QSO sightlines. Sec- 
ond, a skewed frequency distribution of Mg II absorbers is 
also expected along GRB sightlines. None is confirmed in 
empirical data (e.g. Pontzen et al. 2007). We refer the readers 
to Porciani et al. (2007) for a detailed discussion. 

Late-time HST images of the fields around known GRBs 
offer a detailed view of the galaxy environment around known 
foreground Mg II absorbers along the sightlines toward the 
GRBs (e.g. Pollack et al. 2008). Six of the 15 GRBs sight- 
lines in our sample have known strong Mgll absorbers of 
W(2796) > 1 A in the foreground (e.g. Prochter et al. 2006), 
three of which have deep HST images available. Including 
the field around GRB 060418 (Pollack et al. 2008), we find 
that in all four cases the HST images have uncovered at least 
one galaxy at angular separation A 9 < 1" from the afterglow 
sightlines. In contrast, such galaxies at small A 8 are clearly 
absent in fields with no known strong Mg II absorbers (Figure 
22). A summary of the known Mg II absorbers and possible 
candidate absorbing galaxies is presented in Table 4. 

Although only one of the candidate galaxies in Table 4 
has been spectroscopically confirmed (the Mg II absorber at 
z = 0.84 toward GRB 030429), the finding of additional faint 
galaxies along the GRB sightlines strongly disfavors the ab- 
sorbers being intrinsic to the GRBs. In addition, these candi- 
date galaxies exhibit a broad range of optical-infrared colors 
(c.f. the blue colors of galaxies toward GRB 050820 shown in 
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Table 2 and Figure 9 and those found toward GRB 060418 in 
Pollack et al. 2008). Dust extinction does not appear to be a 
uniform factor across different GRB sightlines that results in 
the higher incidence of strong Mgll absorbers. Comparisons 
of the absorbing galaxy properties found along afterglow and 
QSO sightlines (O'Meara et al. 2006; Becker et al. in prepa- 
ration) should provide further insights for understanding the 
differential incidences of Mg II absorbers. The morphology 
of additional faint emission near the absorbing galaxy toward 
GRB 030429 (Figure 6) is suggestive of a lensed event. Ad- 
ditional optical imaging and spectroscopic data are necessary 
to test the lensing hypothesis for this source. 

A conclusive answer to the observed overdensity of fore- 
ground Mg II absorbers along GRB sightlines requires a larger 
sample of imaging and spectroscopic data of the candidate ab- 
sorbing galaxies. Based on the current finding, however, we 
caution that the presence of intervening galaxies at small an- 
gular distances to the GRBs introduces non-negligible con- 
tamination for identifying GRB host galaxies based on imag- 
ing data alone. High spatial resolution images are crucial for 
resolving the host galaxies from foreground absorbers. While 
the properties of the Mgll absorbers are beyond the scope 
of this paper, we emphasize that in the absence of spectro- 
scopic observations it will be necessary to take into account 
line-of-sight absorption-line properties in afterglow spectra 
when evaluating the uncertainty of an imaging identification 
of GRB host galaxies. 

6. SUMMARY 

We present a study of faint galaxies uncovered along GRB 
lines of sight, based on an optical and near-infrared imaging 
survey of the fields around 15 GRBs at z > 2. The GRBs 
are selected with available early-time afterglow spectra in or- 
der to compare ISM absorption-line properties with stellar 
properties. The redshifts of the GRBs span a range from 
Z = 2.04 to z = 4.05. The neutral hydrogen column densities 
of the GRB host ISM span a range from log N(U I) = 16.9 to 
log 7V(H I) = 22.6. Our analysis differs from previous studies 
in that we have obtained a uniform set of photometric data 
from our own imaging survey, reducing systematic uncertain- 
ties in photometric measurements of the host galaxies. 

In addition to the five previously studied GRB host galax- 
ies, we consider new detections for the host galaxies of 
GRB 050820 and GRB 060206. We also place 2-a upper lim- 
its for the rest-frame luminosities of the remaining eight GRB 
host galaxies based on the depths in available optical and near- 
infrared images. Combining early-time, high-resolution after- 
glow spectra and late-time imaging survey of the GRB fields 
allows us to address a number of issues regarding both the na- 
ture of GRB progenitor environment and star-forming physics 
in distant starburst galaxies. The line-of-sight properties un- 
covered in the afterglow spectra have also proven to be valu- 
able for filtering potential contaminations due to foreground 
galaxies. The results of our study are summarized as the fol- 
lowing: 

1 . GRB host galaxies exhibit a broad range of rest-frame 
U V absolute magnitudes spanning from Mab(U) - 5 log h > 
-15 to M AB (UV)-5 log h = -19.8. The distribution of rest- 
frame UV luminosities shows that the GRB host galaxy pop- 
ulation is best described by a UV luminosity weighted random 
galaxy population with a median luminosity of (L(UV)) = 
0.1 L*. Models that include no luminosity weighting are ruled 
out at > 99% confidence level. This result demonstrates 
that GRB host galaxies are representative of unobscured star- 



forming galaxies at z > 2. 

2. There exists a relatively significant correlation between 
UV luminosity and Si II A 1526 transition in GRB host galax- 
ies. A generalized Kendall test including upper limits indi- 
cates that the probability of a positive correlation between 
M uv and W(1526) is nearly 95%. Attributing the observed 
large line widths of W(1526) > 1.5 A to gravitational motions 
of gaseous clouds would require massive halos that are also 
more rare at z > 2. Adopting M uv as a measure of on-going 
SFR and the lack of correlation between SFR and total stellar 
mass in star-forming galaxies at z ~ 2, we therefore interpret 
the observed Muv versus W(1526) correlation as indicating 
that the velocity field observed in GRB host galaxies is driven 
by galactic outflows. 

3. GRB host ISM exhibit a broad range of chemical enrich- 
ment, from < 1/100 solar to ~ 1 /2 solar. No apparent metal- 
licty cutoff is seen in the high-redshift host galaxy popula- 
tion. Similar to nearby galaxies, a tentative trend of declining 
ISM metallicity toward fainter (< 0. 1 L») luminosities is seen 
in the star-forming galaxy population at z = 2-4. The slope 
is steeper than what is expected in some numerical simula- 
tions that incorporate supernovae feedback. Together with an 
absence of molecular gas and the presence of large amounts 
atomic gas, the observed luminosity-metallcity relation may 
be explained by low star formation efficiency in dwarf galax- 
ies. 

4. We measure the interstellar radiation field using resolved 
rest-frame UV morphologies in available HST images. We 
find that the U V radiation field in GRB host ISM spans a range 
over « 35- 350 x higher than the Galactic mean value. The 
strong ISM radiation field observed in GRB host galaxies is 
expected to increase the formation threshold of molecules and 
suppress subsequent star formation, supporting the hypothesis 
that star formation efficiency in dwarf galaxies is reduced due 
to the strong radiation field of existing H II regions. 

5. We examine the galaxy environment in available HST 
images of strong Mg II absorbers found along GRB sightlines. 
In all fields with known strong absorbers, we identify at least 
one faint galaxy at < 1" from the afterglow position. In con- 
trast, such galaxies at small A 9 are clearly absent in fields 
with no known strong Mgll absorbers. The finding of addi- 
tional faint galaxies along the GRB sightlines strongly disfa- 
vors the strong absorbers being intrinsic to the GRBs. Ad- 
ditional imaging and spectroscopic data of the candidate ab- 
sorbing galaxies are necessary to investigate the effect of dust 
and gravitational magnification, but the presence of interven- 
ing galaxies at small angular distances to the GRBs increases 
the ambiguity of identifying GRB host galaxies. It is neces- 
sary to combine high spatial resolution images and early-time 
afterglow spectra for accurate identifications of the host galax- 
ies. 
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TABLE 3 

Summary of known ISM and Stellar Properties of GRB Host Galaxies at z > 2 



Absorption Properties 



Galaxy Properties 









W(1526)° 








M AB (UV) 





M AB (B)f 




rlcltl 


ZGRB 


lOg 1\ yrl 1) 


(A) 


rMn* 
L H J ISM 


[Ml' 
LFeJiSM 


A V 


—5 log h 


(photons/cm /s/Hz) 


—5 log h 


(xlO 9 hr 2 M©) 


GRB000301C 


2.040 


21.2±0.5 










-15.2 ±0.5 








GRB 000926 . 


2.038 


21.3±0.3 


2.40 ±0.14 


-0.17 ±0.15 


+1.32±0.15 


0.15 


-19.63 ±0.07 


3.3 X 10- 6 


-20.26^ 
-19.1 ±0.3 


2.1+40 
— i.y 


GRB 01 1211 . 


2.140 


20.4 ±0.2 


1.35 ±0.19 


>-1.3 






-19.2 ±0.1 




GRB 02 1004 . 


2.329 


19.5 ±0.5 




< -1.0 






-19.83 ±0.07 




-20.38 ±0.15 


-1.9 


GRB 020124 . 


3.198 


21.7 ±0.2 










> —15.1 






GRB 030323 . 


3.372 


21.90±0.07 


0.75 ±0.03 


>-1.26 






-17.6 ±0.1 


8.8 X 10- 7 


18.8 




GRB 030429 . 


2.658 


21.6±0.2 














>-20.1 




GRB 050401 . 


2.899 


22.6±0.3 


2.31 ±0.26 


> -1.57 






>-17.4 




>-19.6 




GRB 050730 . 


3.968 


22.15±0.10 


0.37 ±0.05 


-2.26 ±0.14 


+0.24 ±0.11 


0.00 


> -18.6 








GRB050820A 


2.615 


21.0±0.1 


1.65 ±0.05 


-0.63 ±0.11 


+0.97 ±0.15 


0.08 


-18.50±0.06 


5.6 x 10- 6 


-19.3 ±0.3 


9 +2 5 


GRB 050908 . 


3.343 


17.55 ±0.10 










> -18.9 






GRB 050922C 


2.199 


21.5 ±0.1 


0.52 ±0.05 


-2.03 ±0.15 


+0.60 ±0.10 


0.01 


> -17.1 




> -18.2 




GRB 060206 . 


4.048 


20.85 ±0.10 


1.15 ±0.05 


-0.85 ±0.15 






-17.7 ±0.1 


9.1 x 10- 6 






GRB 060607 . 


3.075 


16.85±0.10 










> -18.3 




> -18.3 




GRB070721B 


3.626 


21.50 ±0.20 










> -19.3 









a Rest-frame absorption equivalent width of Si II A 1526 in the host ISM. The measurements of all but the one for GRB 060206 are adopted from Prochaska et al. (2008). The measurement for 
GRB 060206 is determined using public spectra in Subaru Science Data Archive (Aoki et al. 2006). 

b All measurements but the one for GRB 000926 are based on the observed S II absorption strength. The measurement for GRB 000926 is based on the observed Si II absorption strength, 

because S II transitions are not covered in available moderate resolution spectra of the afterglow. All lower limits are based on the observed Zn abundance in low-to-moderate resolution 

afterglow spectra. The upper limit for GRB 021004 marks the maximum Si abundance prior to an ionization fraction correction. 

c The relative abundances are measured based on sulfur wherever sulfur is available. For GRB 00096, the reported value is based on silicon. 

d The values are derived, assuming the SMC extinction law (Gordon et al. 2003). 

e ISM far UV radiation field estimated from resolved host galaxy images in HST data. 

f The rest-frame S-band absolute magnitude is derived based on available #-band photometry in the observed frame. For sources at z > 3 with detected emission, we infer Mab(B) from the 
measured Mab(UV) and the mean color of starburst galaxies at z = 2-3 from Shapeley et al. (2005). 



n 

cr 

CD 
CD 



24 



TABLE 4 

Summary of Candidate Galaxies Associated with Foreground W(2796) > 1 A Mg II Absorbers 



Mg II Properties Galaxy Properties 



Field z GRB z Mgn W(2796) a AB A8(") Filter 



GRB 021004 2.329 1.38 1.81 > 24.4 Appl. Phys. Lett. 0.3 F606W 

1.60 1.53 >24.4 Appl. Phys. Lett. 0.3 F606W 

GRB030429 2.658 0.84 3.30 20.57±0.05 ss 1.3" H 

GRB050820A 2.615 0.69 2.99 26.30 ±0.05 ss 0.4" F625W 

1.43 1.89 26.20±0.05 a> 1.3" F625W 

GRB060206 3.548 2.26 1.60 26.22 ±0.55 ss 1.0" F814W 

GRB070721B 3.626 3.09 ... 23.7 ±0.1 »* 0.9" 



a Measurements of Mg II absorbers along GRB 021004, GRB 030429, and GRB 050820A are adopted from Prochter et al. (2006). The Mg II absorber toward GRB 060206 is iden- 
tified based on our own analysis of available FOCAS spectra. The strong absorber toward GRB 07072 IB is identified based on a DLA feature by Fynbo and collaborators (private 
communication). 



